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SECTION  1 


INTRODUCTION 

This  report  presents  several  independently  developed  and  tested  algo¬ 
rithms  implemented  in  Pascal  and  FORTRAN  77  for  determining  if  a  point 
falls  within  a  specified  quadrilateral  Area  of  Interest  (AOI).  These  algo¬ 
rithms  represent  several  ways  of  mathematically  solving  the  AOI  point 
inclusion  problem,  and  provide  a  set  of  unclassified  benchmarks  against 
which  existing  system  algorithms  can  be  evaluated.  These  algorithm  analyses 
and  evaluations  are  being  conducted  for  the  Combat  Developers  Support 
Facility  (CDSF)  at  the  United  States  Army  Intelligence  Center  and  School 
(USAICS)  by  the  Algorithm  Analysis  Subgroup  of  the  USAICS  Software  Analysis 
and  Management  Systems  (USAMS)  Development  Task.  It  is  the  ninth  in  a 
series  of  reports  examining  algorithms  used  in  Intelligence  and  Electronic 
Warfare  (I/EW)  systemsr  Each  of  the  other  reports,  listed  in  Appendix  A, 
analyzes  algorithms  from  several  existing'^/EW  systems  that  perform  a 
single  function  and  examine  their  underlying  mathematics.  The  main  func¬ 
tional  areas  studied  have  been  correlation,  geographic  transformations,  and 
direction  finding  and  location  (fix)  estimation. 

In  this  introductory  section,  the  basic  mathematical  and  geometric 
considerations  on  which  the  AOI  algorithms  presented  later  are  based  will 
be  introduced,  and  the  potential  order  dependence  inherent  in  defining  a 
polygon  by  the  location  of  its  vertices  discussed.  The  second  section 
presents  the  algorithms,  and  the  third  section,  an  overview  of  testing 
results  and  algorithm  evaluations. 
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Detailed  test  procedures  and  results  are  given  in  Appendix  C.  The  algo¬ 
rithms  presented  here  were  developed  by  several  different  individuals.  The 
subsection  describing  each  algorithm  is  authored  oy  the  principal  developer 
of  that  algorithm,  with  acknowledgments  to  associates  as  appropriate. 

1.1  UNDERLYING  MATHEMATICAL  CONCEPTS 

There  are  two  basic  mathematical  approaches  to  determining  if  a  point 
is  included  within  the  area  bounded  by  a  closed  polygon.  The  topological 
approach  relies  on  some  notions  of  connectedness  and  the  Wandering  Cow 
Theorem.  The  analytic  approach  uses  the  concept  of  winding  number  and  is 
baaed  on  Cauchy's  Theorem.  These  approaches  will  be  discussed  in  this 
section,  and  algorithms  presented  in  Section  2  based  on  both  methods. 

Consider  the  polygon  P  formed  by  "connecting  the  dots"  in  Figure 
1-1(a)  in  the  order  specified,  giving  Figure  1-1(b).  P  is  a  closed  curve 
(since  the  last  and  first  points  are  the  same)  and  the  region  R  it  end  .js 
is  convex,  that  is,  any  two  points  (e.g.  points  a  and  b)  within  it  can  be 
joined  by  a  straight  line  which  also  lies  within  the  curve.  Also,  the  line 
connecting  any  point  inside  the  closed  curve  P  with  a  point  outside  it 
(e.g.  points  a  and  c)  must  cross  P  exactly  once.  This  is  the  basic  idea 
behind  the  topological  approach  to  determining  if  a  point  lies  inside  or 
outside  the  closed  polygon. 

There  are  formally  two  ways  of  viewing  this  intuitive  result.  The 
polygon  P  is  the  boundary  of  the  region  R,  and  the  set  P+R  is  closed  ("+" 


being  used  to  represent  the  union  operator).  Then  any  connected  set  (i.e. 
any  line  segment  is  a  connected  set)  that  does  not  cross  the  boundary  P  of 
R  lies  either  in  R  or  outside  P+R  (Moore,  Theorem  33).  Or,  viewed  another 
way,  given  the  point  a  in  R  and  c  outside  P+R,  then  P  cuts  between  a  and  c 
(Hocking  and  Young,  Theorem  3-6).  This  means  P  intersects  every  path 
(actually,  every  closed,  connected  subset)  containing  both  a  and  c.  Thus 
the  cow  can  not  wander  from  the  pasture  onto  the  road  without  crossing  the 
fence. 

Neither  theorem  cited  above  requires  R  to  be  convex.  These  theorems 
do  place  restrictions  on  the  topological  space  containing  P  and  R,  but  in 
this  case  P  and  R  are  so  embedded  in  the  Euclidean  plane  that  the  restric¬ 
tions  are  easily  satisfied.  Thus,  at  first  glance,  these  results  would 
seem  to  apply  to  any  polygon.  However,  there  is  a  vacuous  case,  the 
straight  line  (Figure  1-1  (c))  where  they  do  not  help,  because  R  is  empty  so 
contains  no  point  a  and  certainly  can  not  contain  a  connected  set.  But  in 
all  other  cases  the  theorems  apply,  although  for  polygons  such  as  that 
illustrated  in  Figure  1-1(d),  implementing  this  approach  becomes  tedious 
and  tricky. 

The  analytic  approach  is  based  on  Cauchy’s  Theorem  which  says  that  the 
value  of  the  integral  of  an  analytic  function  f(z)  (that  is,  a 
differentiable  complex  function)  around  a  regular  closed  curve  is  zero. 
There  are  some  restrictions  on  the  polygon  P  hidden  in  the  adjectives 
"analytic"  and  "regular."  Certainly  it  is  closed  by  definition,  since  the 
last  point  is  connected  back  to  the  first  point.  The  other  conditions  are 


certainly  satisfied  by  any  convex  polygon  with  three  or  more  vertices,  such 
as  that  developed  in  Figure  1-1(a  and  b).  They  are  also  satisfied  if  the 
polygon  is  merely  star-like  (Figure  1-2).  In  fact,  any  finite  number  of 
smooth  curves  joined  at  their  endpoints  and  enclosing  an  open  region  will 
work  (Figure  1-3).  However,  polygons  such  as  those  illustrated  in  Figure 
1-4  may  cause  trouble  due  to  segments  which  enclose  no  area.  Carefully 
handled,  (a)  can  be  viewed  as  two  disjoint  closed  curves,  for  which  the 
theorem  holds,  but  (b)  can  not. 

Now  using  Cauchy's  Theorem  it  can  be  shown  that,  for  a  point  a  not  on 
the  closed  curve  P,  the  value  of  the  integral  around  P  of  dz/(2-a)  is  a 
multiple  of  (360*^)i.  This  can  be  seen  (but  not  proved)  by  observing  that 
dz/(z-a)  =  d  log(z-a) 

=  d  log!z-a,'  +  i[d  arg(z-a)]; 

the  integral  around  a  closed  curve  of  the  first  term  is  zero  and  arg(z-a) 
increases  or  decreases  by  a  multiple  of  3600. 

Define  the  winding  number  n(a,P)  as  this  integral  divided  by  (360®)i.  Then 
n  is  positive  if  P  is  traversed  counterclockwise,  and 

n(a,-P)  =  -n(a,P). 

Since  R+P  is  a  bounded,  connected,  closed  set, 

n(a,P)  r  0 

for  all  a  not  in  R+P.  Furthermore, 


n(a,P)  =  n(b,P) 


for  all  a  and  b  in  R.  The  only  place  the  winding  number  is  not  defined  is 
on  P  itself.  This  once  again  excludes  Figures  1-1 (c)  and  1-4(b)  from 
consideration.  However,  it  also  means  an  algorithm  based  on  this  principle 
must  check  if  the  point  lies  on  P  for  anv  polygon. 

One  final  result  is  useful  when  designing  AOI  point  inclusion  algo¬ 
rithms.  It  pertains  to  cases  where  n  =  1. 

Let  a  and  b  be  two  points  on  the  polygon,  and  let  P^^jj  be  the  arc  from 
a  to  b,  Pba.  from  b  to  a  (Figure  1-5).  Suppose  that  a  lies  in 

the  lower  half  plane  and  b  in  the  upper  half  plane.  If  P^jj  does  not 
meet  the  negative  real  axis,  and  P^a  ®®®^  ^•^®  Positive  real 

axis,  then  n(0,P)  =  1. 

Since  the  winding  number  is  translation  and  rotation  invariant,  the  condi¬ 
tions  on  a  and  b  are  satisfied  without  loss  of  generality.  This  gives  a 
simple  criteria  which  covers  the  n  equals  0  and  1  cases. 


1.2  ORDER  DEPENDENCE  IN  SPECIFYING  VERTICES 

Throughout  the  above,  the  vertices  were  assumed  specified  in  the  order 
in  which  they  were  to  be  connected.  Can  this  restriction  be  relaxed? 
Unfortunately  not. 


Given  four  points  in  the  plane,  if  one  of  them  is  in  the  interior  of 
the  convex  hull  of  (that  is,  the  smallest  convex  set  containing  the  other 


three  (Figure  1-6  (a)),  the  four  points  do  not  define  a  unique  quadri¬ 
lateral.  In  general,  they  will  define  three  different  ones  (Figure  1-6(b), 
(c),  and  (d))  having  only  the  line  segments  AD,  BD,  and  CD  common  to  all 
three.  Thus  any  point  in  the  convex  hull  that  does  not  fall  on  one  of 
these  three  line  segments  (the  shaded  area  in  1-6(e))  is  or  is  not  in  "the" 
quadrilateral  depending  on  which  quadrilateral  is  chosen. 

Therefore,  to  be  able  to  determine  if  a  fifth  point  falls  within  a 
quadrilateral  defined  by  four  given  corner  points,  either  the  order  (clock¬ 
wise  or  counterclockwise)  of  the  points  must  be  specified,  or  the  resulting 
quadrilateral  must  be  convex  and  the  boundary  of  the  convex  hull  taken  to 
be  the  quadrilateral.  Any  polygon  with  four  or  more  vertices  can  be  made 
concave  by  reordering  the  defining  points  (see  Figure  1-6(f));  remember 
that  these  self-intersecting  polygons  are  mathematically  the  hardest  to 
deal  with.  To  avoid  constraining  the  problem  by  considering  only  convex 
quadrilaterals,  and  it  is  not  always  obvious  that  a  figure  is  not  convex, 
nor  that  easy  to  check  for  convexity,  the  order  of  the  vertices  must  be 


specified. 


Figure  1-5.  An  n(0,P)  =  1  Example 


SECTION  2 


ALGORITHMS 


Before  presenting  the  algorithms  themselves,  a  formal  statement  of  the 
problem  is  needed. 


Statement  of  Problem: 


Given  an  ordered  set  of  n  points  Pi(xi,yi),  P2  ( *2  *  ^2  ^  » 
Pn(Xn,yn)  where  n  may  be  any  finite  positive  integer,  determine  if  a 
specified  point  P(x,y)  is  enclosed  by  the  polygon  formed  by  connecting 
the  points  in  the  given  order  by  straight  line  segments. 

This  polygon  may  be  concave  or  convex  and  may  loop  in  any  fashion. 
The  point  P  may  be  multiply  enclosed  by  the  polygon,  in  which  case  some  of 
the  algorithms  will  determine  the  number  of  enclosures. 


2.1  A  CONVEXITY  ALGORITHM  (QUADTEST) 
Author:  E.  Drell 

Aclcnowledginenta:  R.  Henderson,  CSUN 

E.  Schmidt,  JPL 

Principle  of  the  Algorithm: 


This  algorithm  was  designed  specifically  for  quadrilaterals.  It  is 
based  on  three  observations. 


(1)  Any  path  connecting  a  point  inside  and  a  point  outside  a  closed 
curve  intersects  the  curve  (the  Wandering  Cow  Theorem  discussed 
above ) . 

(2)  Any  quadrilateral  has  at  most  one  point  of  concavity. 

(3)  Any  quadrilateral  (Figure  2- 1(a))  can  be  divided  into  two  trian¬ 
gles  (Figure  2-1(b)). 

The  last  point  allows  the  first  to  be  applied.  Although  in  general  it 
would  be  hard  to  establish  that  any  path  did  or  did  not  cross  the  boundary, 
in  the  case  of  a  triangle  this  reduces  to  showing  that  the  lines  passing 
through  AP,  BP,  and  CP  intersect  BC,  CA,  and  AB,  respectively  (see  proof  in 


Note:  This  algorithm  produces  valid  results  for  test  cases  C-2,  C-5,  and 
C-7 ,  and  is  not  applicable  to  the  other  test  cases. 


I 

f 
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Appendix  D).  Thus  the  following  algorithm  will  divide  the  quadrilateral 
into  two  triangles,  then  use  this  test  to  determine  if  the  point  lies 
within  either  triangle. 


The  Algorithm:  (implemented  in  QUADTEST  in  Appendix  C) 

1.  Test  for  vertices  and  output  message  and  stop  if  point  is  on  a 
vertex. 

2.  If  not,  determine  the  concavity  of  the  figure. 

To  do  this  determine  if  either  the  second  or  forth 
points  of  the  figure  lie  within  the  triangle  defined  by  the  other  3 
points.  Recall  this  can  only  happen  for  one  point.  If  this  condition 
is  true,  then  assign  the  second  vertex  as  the  point  of  inflection  or 
concavity.  If  the  condition  is  false  then  assign  the  first  vertex  as 
the  point  of  inflection.  The  concavity  only  determines  the  way  the 
quadrilateral  is  divided. 

3-  Break  the  quadrilateral  into  two  triangles  based  on  the  inflec¬ 
tion  point  (Figure  2-1  (b)). 

The  default  is  taking  the  first  and  third  vertices 
as  the  center  line.  If  so  indicated  by  the  inflection 


point,  the  second  and  fourth  are  used  instead. 


Test  the  two  triangles  to  see  if  the  point  is  interior  to  either. 


A.  For  each  side  construct  a  line  from  the  vertex  opposite  the 
side  to  the  point  in  question. 

B.  If  this  line  intersects  the  triangle  side  within  the  bounds 
of  the  two  endpoints  of  the  side,  then  the  test  is  successful. 
For  any  point  to  be  interior  to  the  triangle,  each  vertex/point 
line  must  intersect  its  respective  edge.  If  less  than  three 
intersections  are  determined  then  the  point  is  outside  the  trian¬ 
gle.  To  test  if  a  vertex/point  line  and  triangle  side  intersect, 
subtract  the  equation  of  the  line  defining  the  triangle  side  from 
the  equation  of  the  line  constructed  in  (A).  Evaluate  the  re¬ 
sulting  equation  at  the  two  endpoints  of  the  side  yielding  two 
values.  For  a  solution  equal  to  zero  to  exist  within  the  inter¬ 
val  between  the  two  points,  the  value  on  one  side  of  the  interval 
must  be  greater  than  zero  while  on  the  other  side  less  than  zero 
(by  the  Intermediate  Value  Theorem).  Thus  by  multiplying  th^  tvo 
values  together,  if  a  number  greater  than  zero  results,  the 
solution  does  not  fall  within  the  interval  forcing  a  failure  of 


the  triangle  test. 


2^  AN  ANGLE-SUMMING  ALGORITHM  (R_J_G) 

Author;  R.J.  Gardner 
Acknowledgment:  F.  Lesh 

Principle  of  the  Algorithm 

Use  the  point  in  question,  P(x,y),  as  the  vertex  of  all  angles.  Cal¬ 
culate  the  angles  P1PP2,  P2PP3....Pn-1  PPn»  PnPPi  Pn  =  Pt 

Call  these  angles  A^  where  i=1,2,3...n  (i.e.,  A^  is  PiPP2)-  The  order 
and  sign  of  the  angles  must  be  carefully  preserved,  and  the  value  of 
each  angle  must  be  properly  determined  in  the  interval  -180°^  Aj^^l80°« 
If  the  algebraic  sum  of  these  n  angles  (the  winding  number)  is  zero, 
the  point  in  question  Is  not  enclosed  by  this  polygon  as  defined  by 
the  ordering  of  the  points.  If  the  absolute  value  of  this  sum  is 
360^n,  then  the  point  has  been  enclosed  n  times  by  the  polygon. 

The  Algorithm:  (implemented  in  R_J_G  in  Appendix  C) 

Step  1.  Read  in  the  polygon  points  Pi(xi,yi),  P2 (X2 , ys ) .  •  • • » 
Pp(Xn,yn)»  ”  P  resolution  to  consider  two  points  coincident). 

Step  2.  Read  in  the  test  point. 

Step  3-  Initialize  sum  S  to  zero. 

S  =  0 
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Step  4.  For  i  equals  1  through  n  do  the  following;  Translate  the 
origin  of  the  coordinates  to  the  test  point,  P(x,y). 

Xi  =  Xi  -X, 

Yi  =  Yi  -y- 

NOTE:  P(X,y)  =  (0,0). 

Calculate ; 

Mi  =  sqrtCXi+Y^. 

If  Mi  ^  R  output  message  "Point  of  interest  is  on  the  vertex." 
Go  to  Step  8.  This  branch  is  used  to  prevent  division  overflow. 
R  was  chosen  as  0.0001  for  the  purpose  of  this  test  (1  km  in 
10,000  km). 

Step  5.  Do  this  step  for  each  of  the  i  equal  to  1  through  n  with 
n  +  1  =  1. 

Ai  =  C3ign(Xiyi.,.1-Xi.,.iyi)]arccos[(XiXi+i+yiyi+i)/(MiMi^l)] 

As  Ai  approaches  ±180°,  the  testpoint  approaches  the  side.  For  the 
purpose  of  this  test,  ±179.88°  was  used  to  ensure  the  capture  of 
boundary  cases  in  spite  of  round  off  errors,  (i.e.  If  jAi'  >  179.88 
then  output  the  message,  "The  point  is  ON  THE  LINE.") 

S  =  S+Ai 

Step  6. 

E  =  Round[abs(S/360° )]  to  the  nearest  integer. 


2.3  THE  QUADRANGLE  METHOD  (TEST_IT) 

Author:  F.  Lesh 

Principle  of  the  Algorithm; 

Trigonometry  may  be  thought  of  as  the  study  of  the  mapping  of  spatial 
points  in  two  dimensions  onto  a  unit  circle  at  the  origin  of  the  coordinate 
system.  The  angles  are  measured  by  their  arc  length  on  the  unit  circle 
(see  Figure  2-2  (a)),  and  are  positive  in  the  counter-clockwise  sense. 
This  quadrilateral  algorithm  is  similar  to  Section  2.2,  R_J_G,  which  is 
based  on  the  unit  circle.  However,  rather  than  angles,  the  points  are 
mapped  onto  a  unit  square  defined  by  (1,1),  (-1,1),  (-1,-1),  (1,-1)  in  that 
order.  Analogous  to  arc  lengths  as  angles,  call  these  displacements  along 
the  unit  square  "quadrangles."  The  measure  of  these  "quadrangles"  is  the 
distance  along  the  perimeter  between  the  ordered  mapped  points  starting  at 
point  (1,1).  In  Figure  2-2  (b)  the  values  of  the  quadrangles  at  critical 
points  are  given  as  the  eight  values  for  S.  These  critical  points  cor¬ 
respond  to  multiples  of  Tr/4  on  the  unit  circle.  The  origin  of  the  coordi¬ 
nate  system  is  translated  to  the  point  in  question  and  the  zero  quadrangle 
is  defined  by  the  point  (1,1)  Figure  2-3  (b).  Starting  with  the  first 
point  of  the  AOI,  calculate  the  quadrangle  between  consecutive  points 
through  the  last  corner  on  to  the  first  corner.  The  algebraic  sum  (Figure 
2-2  (c))  of  these  quadrangles  divided  by  8  gives  the  number  of  times  the 
point  has  been  enclosed  by  the  polygon.  In  an  analogous  manner,  the  sums 
of  the  angles  divided  by  2Tr  gives  the  number  of  enclosures  of  the  test 


polygon.  This  algorithm  has  the  advantage  over  angle  summing  with  trigono¬ 


metry  in  that  it  uses  no  transcendental  functions,  uses  less  memory  to 
execute,  and  executes  more  rapidly. 

The  Algorithm:  (implemented  in  TEST_IT  in  Appendix  C) 

Given  n  points  on  a  plane,  (Xo,yQ),  (x^  ,y^ ),..., (x^. ,  .y^.i ) ,  and  a  test 
point  (h,k);  determine  if  the  teat  point  lies  inside,  outside,  or  on  the 
boundary  of  the  area  bounded  by  the  lines  connecting  (xQ,yQ)  to  (xi,y-|)  to 

(X2,y2) . (Xn-1.yn-l)  to  (xQ.yo)-  Note  that  Xn=xo,  and  yn=yo- 

Two  functions  must  be  defined. 

S(u,v)=  the  distance  around  the  unit  square  (Figure  2-2  (b)  from  (1,1) 
counter  clockwise. 


If 

Ivliju! 

AND  v>0,  THEN  S(u,v) 

=  1-(u/v). 

If 

lulllvl 

AND  u<0,  THEN  S(u,v) 

=  3+(v/u). 

If 

Ivj^ul 

AND  v<0,  THEN  S(u,v) 

=  5-(u/v). 

If 

1 U  {21 V  1 

AND  u>0,  THEN  S(u,v) 

=  7+(v/u). 

R(w)=  a  function  used  to  calculate  the  difference  between  two  points  as 
seen  from  the  test  point.  A  small  variable  e,  is  introduced  to  ensure  the 
boundary  conditions  against  round  off  errors.* 


•Note;  An  arbitrary  value  for  <  of  0.0001  was  chosen  for  the  "Test  and 
Input"  in  Appendix  C. 
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If  -4+€<W<4-€  R(w)  =  w. 


If  w<-4-€  R(w)  =  w+8. 

If  w>4+€  R(w)  =  w-8. 


Otherwise  test  point  is  on  a  boundary,  and  go  to  next  case. 


Test  points  on  boundaries  (corners  and  sides)  are  considered  to  be 
inside  the  polygon. 


Proceed  as  follows:  first  calculate 


uq  =  XQ-h,  vq  =  yo-k. 

So  =  SCuq.vq)- 
Z  =  0. 


Then,  for  each  i  from  1  to  n  perform  the  following  sequence  of  steps: 


(1) 

“i 

=  Xi-h.  Vi  =  yi-k 

(2) 

If 

((Ui  s  0)  AND  L(Vj 

(3) 

Si 

=  S(Ui,Vi). 

(4) 

°i 

=  R(Si-So). 

(5) 

So 

=  Si- 

(6) 

z  = 

:  Z+Dj^. 

Finally  Z/8  =  the  number  of  times  the  test  point  has  been  enclosed. 
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2.4  HANDEDNESS  ALGORITHM 
Author;  R.  Gardner 
Acknowledgment:  F.  Lesh 

Principle  of  Algorithm: 

To  determine  if  a  point  is  interior  to  a  convex  polygon,  progress 
along  the  sides  and  ask  for  each  side  (pair  of  points);  "is  the  point  in 
question  on  the  left  or  right  of  the  extended  line  segment  of  these  two 
points?"  If  the  point  is  always  on  the  same  side  of  each  segment  for  a 
complete  circuit  around  the  convex  polygon,  the  point  is  interior  to  the 
polygon. 

To  determine  if  a  point  P(x,y)  is  to  the  left  (or  right)  of  a  directed 
segment  Pi(xi,yi)  and  evaluate  the  determinant; 

yi  1 

*2  ^2  ^  ^  ^  (handedness) 

X  y  1 

If  the  determinant  is  0  the  points  are  colinear.  If  h  is  >0  the  P  is 
a  left  turn.  If  h  <0  the  turn  is  right.  See  Appendix  D  for  proofs. 

This  does  not  work  for  a  concave  polygon.  However  in  this  case  one 
can  consider  the  quadrilateral  as  consisting  of  two  overlapping  triangles 
as  illustrated,  (Figure  2-3).  The  first  triangle  (b)  is  formed  by  exclud¬ 
ing  the  concave  vertex.  If  the  h  test  shows  the  point  is  definitely 
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outside  this  triangle  (h  of  different  signs  and  none  zero)  then  the  point 
is  outside  the  AOI.  However,  in  any  other  case  the  second  triangle  still 
must  be  checked.  The  second  triangle  c  is  formed  by  the  concave  point  and 
the  two  points  adjacent  to  it  in  the  ordered  sequence.  If  the  point  is 
definitely  inside  this  triangle  (h^=o  and  of  the  same  sign),  the  point  is 
definitely  outside  the  AOI.  For  the  boundary  cases,  if  both  tests  had  one 
h  =  0,  the  point  is  outside  the  area  of  interest.  If  one  h^  of  one  of  the 
two  triangles  was  zero,  then  the  unknown  point  is  on  the  boundary  of  the 
AOI. 


It  should  be  noted  at  this  point  that  the  handedness  check  is  used  to 
determine  the  nature  of  a  quadrilateral  (see  Figure  2-4  for  the  quadri¬ 
lateral  case).  The  given  h^  are  hi((Xi_i,  yi.i),  (Xi,  y^),  (Xi^i,  yi+i)). 
For  a  convex  quadrilateral  the  hj^  will  all  have  the  same  sign.  For  a 
concave  quadrilateral  one  sign  will  be  different.  For  the  degenerate  case 
of  a  simple  triangle,  one  h  is  zero  and  the  three  remaining  h  are  of  the 
same  sign.  For  the  triangle  with  a  ray  there  is  one  zero  h  and  one  which 
differs  in  sign  with  the  remaining  pair.  The  "bow  tie"  has  two  adjacent  h^ 
greater  than  zero  and  two  less  than  zero. 

The  Algorithm: 

Write  a  routine; 

HAND(x^  ,y^  ,X2,y2.X3,y3,h), 

which  determines  h  for  the  three  points  in  that  order. 
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step  1.  Input  Pi(xi,yi),  P2(X2fy2)»  P3(x3,y3),  and  P4(x4,yi()- 
Step  2. 

Calculate: 

sign  h(P4,Pi,P2)  =  >^1 
sign  h(Pi,P2,P3)  =  h2 
sign  h(P2,P3,P4)  =  h3 
sign  h(P3,P4,PT)  =  h4 

Where ; 

sign  h  =  +1  if  h  >  0, 
sign  h  =  0  if  h  =  0, 
sign  h  =  -1  if  h  <  0. 

This  set  of  {hi^  is  used  to  determine  the  type  of  quadrilateral. 

Step  3.  Input  point  P(x,y). 

Step  4.  If  all  of  the  {h^^}  are  of  the  same  sign  and  non  zero,  it  is  a 
convex  quadrilateral. 

Calculate: 

sign  h(P^,P2,P)  =  Hi 
sign  h(P2,P3,P)  =  H2 
.  sign  h(P3,P4,p)  =  H3 
sign  h(P4,Pi,P)  =  H4 
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If  the  set  {Hj^}  of  the  same  sign  but  non  zero,  output  the  message, 
"Point  is  in  the  AOI!",  and  go  to  Step  8. 

If  the  set  {H^}  ig  of  the  same  sign  but  with  one  zero,  output  the 
message,  "The  point  is  on  the  boundary  of  the  AOI!",  and  go  to  Step  8. 

If  the  set  contains  two  zeros  and  the  other  two  are  of  the 

same  sign  output  the  message  "Point  is  on  a  vertex  of  AOI!",  and  go  to 
Step  8. 

Otherwise  output  the  message  "Point  is  outside  the  AOI!",  and  go  to 
Step  8. 

Step  5.  If  no  hj^  are  zero  and  one  of  the  h,  say  hj,  is  of  a  different 
sign  than  the  other  two,  it  is  a  concave  quadrilateral.  From  the  tri¬ 
angle,  formed  by  excluding  this  Pj^  calculate  the  following; 

sign  h(Pj+i,Pj.^2.P)  =  H'l 
sign  h(Pj^2.Pj+3,P)  =  H'2 
sign  h(Pj^3,Pj+1,P)  =  H*3 

If  any  two  of  these  are  different  in  sign  (excluding  a  zero) 
output  the  message  "Point  is  outside  AOI!",  and  go  to  Step  8. 

If  two  of  the  are  zero  output  the  message,  "Point  is  on  a  vertex 
of  AOI!",  and  go  to  Step  8. 

Otherwise,  for  the  triangle  formed  by  the  Pj  (concaved  point)  and  its 
adjacent  points  (Pj_i  and  Pj+i^  calculate; 

sign  h(Pj_^,  P)  =  H"i 

sign  h(Pj,  Pj^.1,  P)  =  H"2 
sign  h(Pj^,,  Pj_^,  P)  =  H"3 
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If  the  are  all  of  the  same  sign  and  not  zero,  or  if  both  {H'l  and 
{H"}  contain  a  zero,  output  the  message,  "Point  is  exterior  to  t^ 
AOI!",  then  go  to  Step  8. 


If  two  are  zero,  output  the  message,  "The  point  is  on  a  vertex  of 
the  AOII",  and  go  to  Step  8. 

If  either  set  or  {H"^}  contains  a  zero,  but  a  zero  is  not 

contained  in  both,  and  {H"}  contains  no  change  in  sign,  output  the 
message,  "The  point  is  on  the  boundary  of  the  AOI!",  then  go  to 
Step  8. 

Otherwise,  output  the  message,  "The  point  is  inside  the  AOI!",  and  go 
to  Step  8. 

Step  6. 

If  only  one  h  (hj)  ia  zero,  (Figure  2-4(c)  or  (d))  form  a  triangle  by 

excluding  the  point  Pj;  i.e.  Pj+i i Pj4.2» Pj+3 

where  j  +  n  is  cyclic  in  the  range  1,2,3f4.  Read  in  the  P(x,y). 

Then  calculate; 

sign  h(Pj^^,  P)  =  H’l. 

sign  h(Pj,.2,  Pj+3,  P)  =  H'2- 
sign  h(Pj^3,  Pj^,,  P)  =  H'3. 

Where  the  H'j,  are  non  zero  and  of  the  same  sign  output  the  message 
"The  point  ia  in  the  AOI!". 

Where  there  is  one  zero  and  the  other  two  are  of  the  same  sign 
output  the  message  "Point  is  on  the  boundary  of  the  AOI!". 

Where  there  are  two  zero  are  non  zero  and  of  the  same  sign  output 
the  message  "Point  is  on  vertex!",  and  go  to  Step  8. 


Where  none  of  these  conditions  exist  output  the  message,  "Point  is 
outside  the  AOIl". 


Go  to  Step  8. 

Step  7. 

Otherwise  output  the  message,  "Illegal  input!"  and  go  to  next  step  9. 

Step  8.  Output  the  message  "Are  there  more  points  to  Input?".  If 
answer  is  yes  go  to  Step  4. 

Step  9.  Output  the  message,  "Are  there  more  areas  of  interest  so 
specify?".  If  answer  is  yes  go  to  Step  1. 

Step  10,  Exit, 
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(a)  CIRCULAR  ARCS 


(b)  QUADRANGLES 


(a)  CONCAVE  QUADRILATERAL 


(b)  FIRST  TRIANGLE 


(c)  SECOND  TRIANGLE 


Triangulation  of  a  Concave  Quad 


(a)  A  CONVEX  QUADRILATERAL 

4 


(b)  A  CONCAVE  QUADRILATERAL 

4 


(c)  A  SIMPLE  TRIANGLE  (DEGENERATE  QUADRILATERAL) 


(d) 


(e) 


A  TRIANGLE  WITH  A  RAY  (DEGENERATE  QUADRILATERAL) 
4 


(f)  A  LINE  SEGMENT  (DEGENERATE  QUADRILATERAL) 

• - m - >—  - 


Figure  2-4.  All  Possible  Quadrilaterals 


SECTION  3 


TEST  AND  EVALUATION 

The  algorithms  presented  in  Sections  2.1,  2.2,  2.3,  and  2.4  were 
implemented  and  tested  on  a  VAX -11/780  using  Pascal  and  FORTRAN.  The  code, 
test  quadrilaterals  and  points,  test  results,  and  specific  comments  on  the 
performance  of  each  algorithm  are  given  in  Appendix  C.  The  following  is  a 
survey  of  those  findings. 

Two  main  categories  of  performance  emerged  in  designing  and  conducting 
the  tests: 

(1)  Geometric  robustness. 

(2)  Timing  and  efficiency. 

The  first  is  a  descriptive  as  well  as  evaluative  category.  For  example,  in 
some  applications  points  lying  on  the  polygon  itself  may  be  considered 
outside  the  enclosed  region  and  for  other  applications  inside  it.  Table  3- 
1  presents  geometric  robustness  for  the  four  algorithms.  Note  that,  al¬ 
though  in  theory  some  of  the  algorithms  work  for  any  polygons  with  three  or 
more  sides,  only  quadrilaterals  were  used  in  the  testing.  The  three 
classes  tested  were  convex,  concave,  and  degenerate.  Both  design  and 
analysis  of  the  algorithms  indicate  that  for  polygons  with  more  sides  the 
algorithms  presented  in  Sections  2-2  and  2-3  should  work  in  all  but  (possi¬ 
bly)  the  degenerate  cases. 

Questions  of  timing  and  efficiency  arose  primarily  in  five  ways. 
Three  of  these  are  fairly  independent  of  the  specific  application: 
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(1)  How  many  special  cases  had  to  be  checked. 

(2)  How  many  subproblems  were  solved. 

(3)  Whether  trigonometric  functions  were  used. 

Table  3-2  addresses  these.  The  other  major  efficiency  issues  are  quite 
application  dependent.  The  order  in  which  the  various  tests  are  performed 
in  an  algorithm  depend  on  the  characteristics  of  the  most  likely 
quadrilaterals  and  input  points.  Thus  the  following  two  considerations  can 
significantly  affect  efficiency: 

(4)  Whether  to  test  for  wild  inputs. 

(5)  When  to  test  for  degenerate  quadrilaterals. 

Table  3-1.  Geometric  Robustness 


Algorithm  from  Section: 
Name 


2.1  2.2 

QUADTEST 


2.3  2. 

TEST  IT  HANDEDNESS 


o 


Convex  Interiors 

Concave  Interiors 

Boundaries  (Non  Degenerate) 

Degenerate  Quadrilaterals 

- 

Quadrilaterals 

+ 

N-sided  Polygons* 

- 

+  handles 

-  does  not  handle 

*  has  not  been  proved  in  this 

report 

These  questions  are  also  covered  in  Table  3-2.  For  any  applications  the 
trade-off  must  be  made,  whether  there  are  enough  wild  inputs  to  take  time 
testing  for  them,  and  just  how  often  the  quadrilaterals  will  be  degenerate. 
The  table  is  meant  only  as  a  guide  for  applications. 


Table  3-2.  Efficiency  Considerations 


Algorithm  from  Section: 


QUADTEST  R_J_G  TEST_IT  HANDEDNESS 


Subproblems 
Trigometric  Functions 
Wild  Inputs 
Degenerate  Polygons 


Language 


-  does  not  have 


0  comes  late 


♦  comes  early 


F  FORTRAN 


P  PASCAL 


vi-.v 

>•,  *  V ^  w  -I  - 
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APPENDIX  C 


IMPLEMENTATION  AND  TESTING 
by 

Nick  Covella  and  Bruce  Pardo 

For  each  of  the  following  programs  and  test  results,  refer  to  Figures 
C-1  through  C-7.  The  solid  lines  reflect  the  input  quadrilaterals.  The 
X's  indicate  the  locations  of  the  input  test  points.  The  handling  of 
points  on  boundaries  (corners  and  sides)  is  described  in  Section  2. 

There  are  four  coded  programs  (2  in  Pascal  and  2  in  FORTRAN)  that  have 
been  designed  and  tested  by  the  individuals  who  took  part  in  the  design  of 
the  respective  programs.  The  four  programs  are: 


1. 

QOADTEST 

(PASCAL) 

2.1 

2. 

IL.J_G 

(FORTRAN) 

2.2 

3. 

TEST_IT 

(FORTRAN) 

2.3 

4. 

HANDEDNESS 

(PASCAL) 

2.4 

Each  of  these  algorithms  used  different  mathematical  methods  -to 
determine  whether  a  point  lies  inside,  outside,  on  a  corner,  or  on  a 
boundary  of  the  input  figure.  These  methods  are  discussed  in  Section  2. 
The  analyses  of  the  generated  results  appear  below. 


PROGRAM  QUAD:  (SECTION  2.1) 

This  program  used  the  idea  of  triangles.  If  the  line  connecting  the 
test  point  with  each  vertex  intersects  the  line  of  the  opposite  side  of  the 
triangle  at  a  point  within  the  confines  of  the  triangle,  and  this  condition 
holds  true  for  all  three  cases,  then  the  test  point  is  in  the  interior  of 
the  triangle  in  question. 

For  non-degenerate  quadrilaterals,  this  algorithm  performed  well.  It 
does  not  properly  handle  the  bow  tie  or  degenerate  cases  for  which  three 
points  fall  on  the  same  line. 

PROGRAM  RJG:  (SECTION  2.2) 

This  program  used  the  idea  of  angle  summing.  The  point  in  question  is 
used  as  the  reference  point,  then  the  angles  formed  by  this  point  and  the 
points  of  the  figure  are  summed.  Direction  of  counter-clockwise  or  clock¬ 
wise  makes  no  difference  so  long  as  the  points  tested  are  in  sequence.  The 
value  of  the  final  summation  is  either  0  or  a  multiple  of  360.  If  the 
value  is  0,  then  the  point  in  question  is  outside  the  figure,  otherwise  it 
is  inside  the  figure. 


Based  on  the  testing  here,  this  algorithm  seems  very  robust. 


PROGRAM  TEST_IT:  (SECTION  2.3) 


This  algorithm  used  the  idea  of  a  unit  square.  Function  Scale  calcu¬ 
lates  the  distance  around  the  input  figure  starting  at  (0,0)  and  travelling 


counter  clockwise.  Rcalc  calculates  the  angular  difference  between  two 
points  as  seen  from  the  teat  point.  The  result  gives  the  number  of  times 


the  test  point  was  circled.  For  a  closed  quadrilateral  that  had  no  inter¬ 
nal  intersections,  if  the  test  point  was  inside  the  figure  then  the  en¬ 
closed  value  was  1;  otherwise,  it  is  zero. 


Overall,  this  algorithm  seems  very  robust.  It  tested  each  input  point 
correctly  and  quickly  and  could  be  used  in  a  real-time  environment. 


PROGRAM  HANDEDNESS:  (SECTION  2.4) 


This  algorithm's  basic  principle  only  does  not  apply  to  "bow-tie" 
quadrilaterals.  One  should  progress  along  the  perimeter  in  the  direction 
specified  by  the  ordered  points.  With  each  pair  of  the  perimeter  points, 
follow  along  the  direction  specified  by  their  ordering,  connect  two  con¬ 
secutive  points,  form  a  segment,  then  check  to  determine  if  the  test  point 
is  on  the  right  or  left  of  the  segment.  If  the  answer  is  sometimes  "left" 
and  with  other  point  pairs  "right",  the  point  is  exterior  to  the  convex 
quadrilateral.  If  the  answer  is  always  the  same  (left  or  right),  the  point 


is  interior  to  the  convex  polygons.  In  the  case  of  the  concave  quadrila¬ 
teral,  the  convexity  limitation  is  avoided  by  dividing  the  quadrilateral 
into  two  triangles  and  applying  this  principle  to  each  of  the  triangles. 


Though  this  algorithm  appears  very  robust,  the  logic  is  very  compli¬ 
cated  when  handling  the  concave  case.  Thus  it  becomes  impractical  to  im¬ 
plement  in  the  field,  particularly  when  compared  to  the  other  algorithms. 
The  difficulty  of  expanding  this  algorithm  increases  geometrically  with  the 
number  of  sides. 


In  the  following  seven  diagrams  (Figures  C-1  through  C-7),  the  lines 
represent  the  outline  of  the  input  figure,  the  crosses  (X’s)  represent  the 
input  test  points,  and  the  encircled  crosses  are  special  test  points;  the 
vertices  of  the  input  figure. 
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A  4  SIDED  PiJLVGON  GIVEN  THE  4  POINTS  IN  ORDER  AND  THE  COORDINATES 
OF  THE  POINT  IN  QUESTION.  IT  IS  DESIGNED  TO  HANDLE  CASES  OF  BOTH 
CONVEX  AND  CONCAVE  FIGURES. 
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*  1 


TYPE 


POINT  =  ARRAY  LI..  23  OF  REAL; 
TRIANGLE  =  ARRAY  Cl,  .  33  OF  INTEGER; 

VAR 


A 

P 

Tl.  T2 
CENTER 
LFLAG.  I..  J.  K 
INFLECTION 
RESPONSE. REPLY 


ARRAY  C 1.  43  OF  POINT; 

POINT; 

TRIANGl-E; 

INTEGER; 

INTEGER; 

INTEGER; 

CHAR; 
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FUNCTION  TRIANGLETEST  (V1.V2.V3  ;  INTEGER;  P  :  POINT)  :  BOOLEAN; 

(*  THIS  ROUTINE  TESTS  IF  THE  LINE  CONNECTING  THE  POINT  P  WITH  EACH 
VERTEX  INTERSECTS  THE  LINE  OF  OPPOSITE  SIDE  OF  THE  TRIANGLE  AT  A 
POINT  WITHIN  THE  CONFINES  OF  THE  TRIANGLE.  IF  THIS  CONDITION 
HOl.DS  TRUE  FOR  ALL  3  CASES  THEN  THE  POINT  IS  IN  THE  INTERIOR  OF 
THE  TR3ANG1..E  DEFINED. 

*) 


VAR 

VALl.  VAL2..  LX.  RX,.  M 
I , COUNT 


V,  Si  ,  s 

2,  TEMP 

BEGIN 

§ 

Cl  DUN  r 

•  ^  0; 

i 

V  .  = 
Si  = 

ACVl  3, 
ACV23, 

32  :  = 

AlV33; 

FOR  I 

=  1  TO 

r  W 

1 

(  ♦FOP 

THE  T 

uCC 

URS  AT 

REAL; 

INTEGER, 

POINT; 


3  DO  BEGIN 

HREE  SIDES  TEST  IF  THE  INTERSECTION  OF  THE  2  LINES 
A  POINT  BETWEEN  THE  LEFT  AND  RIGHT  ENDPOINTS  OF  THE 
USE  THE  INTFRMErjIATF  VALUE  THEORM  TO  DETERMINE  IF  THE 


SI'iE. 

SOLUTION  EXISTS  IN  THE  INTERVAL 


C-10 


IF  (  31C13  O  S2>;U  ■) 

TMFN  M  :=  (  (  S 1  C  2  ] -S2l2  3  )  /  ( S 1  C  1  3-S2C  1  3  >  )  C^SLQPE-*-) 

EL2E  M  :=  '79999;  (^8LuP2  IS  INFINITY  SINCE  YERTICALi») 

CASE  WHERE  THE  VERTEX  POINT  AND  THE  POINT  IN  QUESTION 
HAVE  THE  SANE  X  VALUE  ONLY  THE  Y  COORDINATE  NEEDS  TO  BE 
TESTED.  A  LINEAR  FORMULA  SIMILIAR  TO  THAT  IN  THE  LINEARITY 
PROCEDURE  13  USED  AND  A  COMPARISON  IS  DONE  USING  THE  TWO 
DIFFERENT  Y  VALUES  WHILE  HOLDING  THE  X  CONSTANT. 


IF  VC13=PC13  OTEST  FOR  SLOPE  OF  INFINITY^) 

THEN  BEGIN 

VALl :  =tvr23-Sl';23  ■;  -  ( M-+ c  VC  1  3-Si  C  1  3  )  ) ; 

VAL  2:  =(DC23 -Si  l23  ;  -  (  N-X  VC  i  3-31  L  1  3  )  ) , 

IF  (ADSLVrtLi)  >  ABSiVALD))  THEN  COUNT  .  =COUNT  +1, 

END  (  *Ti  !£N-* 

EL'OE  BEGIN 

LX:=SiC13;  RX:=S2C13, 

VALi:  =(VC23-32l2j  )-(M^(LX  -32113)  )  ■»■  (  (  ( VC2  3-P  C23  )  /  (  VC  1 3 -P  C  1  3  .'  )'^-(LX-VC 

VAL2;  =(  VC23-31  C2  3  (RX  -SiCl  3  )  )-*-  (  (  (  VC2  3-P  C23  )  /  (  VC  1  3-P  C  1 3  j-XR.i-Vi 

IF  (VAt  1  >  VAL2)  C=  0  THEN  COUNT  ;  =CGUNT  -f  1; 

END  ( <»EL3i-A- ) ; 


.'A 


TEMP  : =  V; 


V 

3 1 

32 


=  SI; 

=  S2i 
-TEMP; 


(♦SWITCH  TO  THE  NEXT  EDGE  OF  TRIANGLE*) 


END  (♦FOR*); 

IF  COUNT  =  3  THEN  TRIANGLETE3T  :=  TRUE  (♦POINT  WITHIN*) 
ELSE  TRIANGLETE3T  :=  FALSE; 

iND. 
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FUNC'l  ION  CONCAVITY 


INTEGER, 
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|♦THI3  ROUTINE  TEST  TO  SEE  IF  THE  FIGURE  IS  CONCAVE  OR  CONVEX  IT 
TESTS  BY  USE  OF  PERMUTAT IONS,  IF  ANY  OF  THE  VERTEX  POINTS  LIES 
WITijIN  THE  TRIANGLE  DEFINED  BY  THE  OTHER  3  POINTS.  IF  CONCAVITY 
13  determined  THEN  THE  FUNCTION  RETURNS  THE  INTERNAL  POINT  AS  THE 
POINT  OF  INri_ECTICN.  IF  THE  FIGURE  IS  DETERMINED  TO  BE  CONVEX 


VER  !  r.  f 


IS  -ETURNED  BY  DEFAULT, 


;  iN 


IF 


END, 


^  (TRIANGLETESTi  1, 2,  3,  AiX'I  )^TRUE)  OR  (  TR  I ANGLETEST  (  3,  4,  1,  AC23  '  =TRL 
THEN  CONCA'-'ITY  ■  =2 
r  L3E  CONCAVITY  =1, 


P70CFDUPE  ENTEP'-'AuE, 

'.♦TAr';  RROCED'.'RE  IS 
^  ♦OF  'VO  I  NTS-*) 


"HE  PRCMP  .  routine  TO  tNTtR  TP-E  VALUbS  INTO  T^t.  aR.R.-'Y^ 

c-n 


U  TLAG  INiiGtH; 

BEGIN 
WRI  I'EI  N; 

WRITE! -N  (  'ENTt-:?.  E/sCH  OF  THE  T  POINTS  IN  ORDER 

WRITE'*. N  <  'ORDER  IS  DETERMINED  BY  THE  CONSECUTIVE  POINTS  BEING  JOINED  BY  A  SIDE 
WRITE'*  N  ■;  'ENTER  E/'-.CH  OF  THE  4  POINTS  AS  AN  ORDERED  PAIR'); 

WR  I  fE;  N; 

F.3R  I:  =  1  TO  4  DO  BEGIN 

WRITE  ('ENTER  T'HF.  X  COORDINATE  OF  VERTEX  '.I;  2.  '  '); 

RC/-.DLN  (AC  I.  11)  . 

WR.TE  ('EMfEK  the  Y  COORDINATE  OF  VERTEX  '>1:2,  ' 

RCAOl.N!  AT  I,  2  i  )  . 

Wt  ;  (  n-ILM; 

KNt); 

WR  I  I  Fi  N; 


r*.)NCT;ON  NEWPi  -  INTtGER- 

('^mi:-;  ROUTINE  ENTERS  THE  POINT  TO  BE  CONSIDERED  AND  DOES  ♦) 

(  <-VKR  *  F I  Cat  [  ON  cHECaING  on  0XI3i'’ING  COORDINATES*) 

BEGIN 

NF'WP  r :  ”0; 

WRiTE  ('ENTER  THE  x  COORDINATE  OF  THE  POINT  P  '), 

KC;'..DLN  (PClD)i 

WRITE  (' enter  the  Y  COORDINATE  OF  THE  POINT  P  '); 

KOADLN  <RCH3); 

WR ( TELN; 

WR (TELN; 

FOh  1:^1  I'O  4  DO 

[F  (  (ACI,  i:=PCn)  .AND  (ATI,  2I=PC23)  )  ( *CHECK  FOR  A  VERTEX  POINT*) 
T'HEN  NEWP  i  :  —  1, 


BEGIN  MAIN  ROUTINE 


aE’j  i  N 

REPEAT  BEGIN 
ON  rFRVALS, 

wr.  .TELN  (  'El'rER  THE  POINT  UNDER  CONSIDERATION  AS  AN  ORDERED  PAIR 
WK i  TELN, 

REr EAT  BEG(N 

LFLAG  -NEWPT. 

IK  LFLAc  i  THEM  BEGIN 
END  (  *END  I  !*  *  ) , 

C.a.SE  LKLAC  OF 

0:  BEGIN  ( *NON  VERTEX  CASE*) 

INFLECTION  :  -'ONCAVITY, 

'ICll:  --INFLECTION,  (*OEFINE  TRIANGLES  BASED  OiN 

TiL2j  ^INFLECT  [ON-*-! ,  (*POINTS 

*'’1131  =  I  i'iFLEC  T 1 0N-t-2i  r-io 


Ti 


T2C13; = INFLECT I ON; 

T2  C  2  3  :  =  I NFLEC  T 1  aN-t-2i 
IF  INFLECTiai'i  ^  2 

THEN  T2C33;  =•  INFLECTION-! 

ELSE  T2i;33:  =INFLECTION  +  2; 

ir  (TRlANGLETESTiTlCU,  T1C23,  T1C33- P>  =  TR  I ANGLETE3T  ( T2C  1 3 ,  TIC  2 
TMFM 

IF  TR  lANGl  ETE3T(T1  C  13,  T1C23,  T1C33,  P)  THEN 

WRITELNI  'THE  POINT  IS  AN  INTERIOR  POINT  OF  THE  GUADRA'-A 
ELSE 

WRITELN(  'THE  POINT  IS  AN  EXTERIOR  POINT  OF  THE  GL'ADRALa 

ELSE 

WRITELNC  'THE  POINT  IS  AN  INTERIOR  POINT  OF  THE  G'JADRALaTER 
END  ^^*0  CONDITION  OF  CASE*); 

I  .  WRITELN  ('THE  POINT  WAS  A  VERTEX  OF  THE  QUADRALATERAL  '  ) ; 

END  (*CASE*)i 

WRITE!  Ni 

WRITE  (-another  point  TO  TEST  WITHIN  THIS  QUADRALATERAL*  Y/N  ■ 
READL.N  I  REPLY), 

WRITELN. 

END  (*REP!_V  REPEAT*); 

UNTIL  REPLY 
WRITELN; 

WRITE  {  'ANOTHER  GUmDRAL  ATERAL  TO  TEST*-  Y/N  '), 

UEADLN  (RESPONSE); 

WRITELN; 

'END  (*fiE3P0NSE  REPEAT*); 

UNTIL  RESPONSE  'O  'Y  '; 

END, 
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CHECKING  THE  INPUT  (X,  Y)  POINTS 


INPUT  X 
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CHECKING  THE  INPUT  (X,  Y)  POINTS 


INPUT  X 
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INPUT  Y 
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0.0000 

5.0000 
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INPUT  Y 
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ACTUAL  RESULT 
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OUT 

OUT 

OUT 

OUT 

OUT 

CORNER 

IN 

OUT 

OUT 

OUT 

CORNER 

CORNER 


L** ' 


- 1.  V- 


QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


C-7  0.0000  0.0000 

10.0000  0.0000 

6.0000  8.0000 

2.0000  4.0000 


CHECKING  THE  INPUT 

(X,  Y)  POINTS 

INPUT  X 

INPUT  Y 

GENERATED  RESULT 

ACTUAL  RESULT 

4.0000 

4.0000 

IN 

IN 

4.0000 

1.0000 

IN 

IN 

2.0000 

3.0000 

IN 

IN 

6.0000 

3.0000 

IN 

IN 

8.0000 

3.0000 

IN 

IN 

5.0000 

0.0000 

IN 

ON  THE  LINE 

2.0000 

7.0000 

OUT 

OUT 

8.0000 

7.0000 

OUT 

OUT 

6.0000 

2.0000 

IN 

IN 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

O.QOQQ 

IN 

ON  THE  LINE 

4.0000 

0.0000 

IN 

ON  THE  LINE 

6.0000 

0.0000 

IN 

ON  THE  LINE 

10.0000 

0.0000 

CORNER 

CORNER 

6.0000 

8.0000 

CORNER 

CORNER 

3.0000 

3.0000 

IN 

IN 

4.0000 

6.0000 

IN 

ON  THE  LINE 

7.0000 

0.0000 

IN 

ON  THE  LINE 

6.0000 

4.0000 

IN 

IN 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

IN 

IN 

12.0000 

0.0000 

OUT 

OUT 

0.0000 

5.0000 

OUT 

OUT 

7.0000 

5.0000 

IN 

IN 

C-20 


'I 


Program  R_J_Q 


This  program  reads  in  -Ix.g)  values  a  three  (Quadrilateral 
having  4  corners.  A  test  point  (x,y)  is  checked  to  deter Tiine  i- 
the  inputted  point  falls  in  the  quadrilaterals 

Algorithm  designed  by:  Robert  J.  Gardner,  J  P  L. 

Implemented  by:  Nick  Covella,  J  P  L 

Dated:  March,  1985 


Real*S  XArr ay < 1 , 4 ) , YAr ray ( 1 , 4 ) 

Real*9  Ne'iiXArray<l,4).  Neu»YArr  ay  <1,45 
Real*S  InXcoord,  InVcoord 
Char ac ter* 1  Resp 1 . Resp2. Flag 
Integer  Quadcntr 

C  *****■<*•■»■*•«•■»■**■«•  *******»****  ♦«■*******•*•»•**■»••«•*•»••«••«•■«■■»■■»■*  **•«••»•■»■*■«■■«■■»••»■  H-*********** 
C* ■**•»•<-*•«■*•»■****><•  *****•«■  •»*«■•«!•♦  MAIN  *«■  PROGRAM  **»-»-***#******  «■*** 


Cpen^Unit  =  1 ,  Fi  1  e- 'R  j  g .  Da  t ' ,  Status  =  'Neiu') 

Open<'Jnit  =  2,  F 1 1  e- '  P  t  s  Da  t ' ,  Status  =  'Old') 

Call  Skip'l) 

Call  Skip<2) 

Call  Pr omp t 1 < Re s p 1 ' 

Quadcntr  =  0 

Do  While  (CRespl  EQ.  'Y')  .OR.  <Respl  EQ.  'y')) 

Quadcntr  =  Quadcntr  +  1 

OpendJnit  =  3,  Fi  lo= 'Test.  Dat '<  Status  *  'Old') 

Cali  Initialiie<XArray, YArray, NewX Array, NewYArray ) 

Call  Load Arrays<X Array. YArray ,  Quadcntr) 

Call  Unl.oa‘‘4rrays(XArray.  YArray.  Quadcntr) 

Call  Pr omp t2 ( R esp 2.  InXcoord,  InYcoord) 

Do  While  ((Resp2  EQ.  'Y')  OR.  (Resp2  EQ  'y')) 

Call  TestPtf InXcoord, InYcoord, XArray, YArray, Flag ) 
Call  Tran5form(  IiiXcoord,  InYcoord,  XArray,  YArray, 
NeuiXArray,  NewYArray  ) 

Call  Te  s  t_I  t  ( NeiiiX  Arr  ay .  NeujY  Arr  a  y  ,  InXcoord,  InYcoord) 
GO  TO  29 

Call  Pr omp t2 ( R 0 s p 2, InXcoord. InYcoord) 

If  ((Re5p2  .  EQ  'N')  OR  (Re5p2  EQ  'n'))  Then 
C 1 o  5  e ( Un i t  =  3 ) 

End  If 
End  Do 

C.dll  Pr  omp  1 1  (  R  e  s  p  1  ) 

End  Do 
Wr  i  to  (  5,  55  > 

C 1 0  se  <  Un 1 1=  1  ) 

C  1  ose ( Un 1 1  =  2 ) 

Format'T2,  'Normal  Completion  oF  Program  by  N  Covella') 

Stop 

End 


C  *************** 


•  J- 


Subroutine  Tesirt'llnXcncrd.  InYcoord.  XArray.  r' Array.  Flag' 
This  Subroutine  checks  if  the  pome  is  one  of  the  corners 


n  o  o  o  o  o  C)  o  o 


Realms  XArr  a  ij  (  1 .  4 )  >  YArr  ay  ( 1  j  4 ) 

Rea  1*3  InXcoord-InYcoord 
Character*!  Riag 

Flag  *  'F' 

Do  to  I  =  1-  4 

[P  vInXcoord  .  EQ.  XArray(l.I))  Then 
If  (InYcoord  .  EQ.  YArray(l- I))  Then 
Flag  =  '1' 

End  If 
End  I  f 

10  Continue 

Return 

End 

C ****»>******** ****** *******  *************** ******************** ********* 
C ********* ******************************************* ****************** 

Subroutine  Initiali  2u(XArray.  YArrayj  Neu.iXArray<  hieixiYArray  ) 

This  subroutine  initializes  the  tujo  arrays  (i.  e  quadrilateral 
and  Analyst)  to  zero  su  that  any  unknown  values  may  be  detected 

Real*S  XArray ( 1 .  4 ) ,  YArray ( 1 / 4 ) 

R  ea  1  *3  NewXArray  <  1 ,  4  ) .  Neu/YArray  (1,4) 

Integer  Index 

Do  10  Index  =  1,4 

Xarray ( 1 , Ind e X )  =  0 
Yarr ay (  1 ,  I nd e X )  -  0 
NewXarr ay ( 1 , Ind e X )  =  0 
NewYarray ( 1 ,  I nd e X  )  =  0 
C  0  n  1 1  n  u  ? 

Return 
End 


Subroutine  LoadArrays(XArray, YArray, Kount) 


This  subroutine  loads  the  3  quadr i lateral s  with  (x.y)  coordinates 
for  ea'ih  corner  of  the  corresponding  quadr  i  lateral  s. 


Rea (*9  XArray ( 1 ,  4 ) ,  YArray ( 1 ,  4 ) 
Intnger  Kount 


W  r  1  t  e  (  5 ..  0 1  ) 

Wr  i  e  (  1 ,  0 1  ) 

Do  to  Index  =  1,4 

Write(5,97)  index, Kount 
Held  (2,99)  XArrjy(i,  Index) 
Ur  1 1  e  (  5 ,  99  ;  I  n  d  e  x  ,  K  o  u n  t 
Head  (2,99)  >'  A  r  r-  j  y  (  1 ,  I  n  d  e  x  ) 
Continue-  p  -j-i 


10 


Oi  Forma  ‘  ■'  ) 

97  Forma?(Tl/  '  Input  X_coordinate  W',  I1/T23/ 

;■  '  for  Quadr  i  latera  1  W',  II) 

9G  Formardl/"'  Input  Y_co ordinate 

:■  '  for  Quadr  i  lateral  #'»  11) 

99  Forma t<F9.  4) 

Return 

End 

C  » -it -i*  ^ -s*  -S- ->■»  rt  -a- -S- ^ -a- -5- -s- -» -s- ■» 

C  *  4^  a- -> ->  4»  ^  41- a- 41- -a- -a- ■» -> ->  ■a' a- ■» -a- T*  ■» -a- 4*.  ij.  i».  s*.  #  4*. »  *  4>.  *  4*.  #  *  » *  4(. »  »  *  » 4*. »  ♦  ■»  *  » -a- *  ■»  4}.  ij.  *  4*.  #  4<. 

Subroutine  UnLoadArraysCXArray  -  YArray  /  QuadNo) 

C 

C  This  subroutine  displays  the  inputted  values  of  the 

C  (x,y)  coordinates  for  each  quadrilateral  and  the 

C  (  .t .  y  )  coordinate  that  is  to  be  checked  luhether  or  not  it 

C  is  in  any  of  the  three  quadrilaterals. 


Realms  XArray(l<  4), YArrayCl.  4) 

Integer  QuadNo,  Index 

UJr  i  te  (  5/  96  ) 

Wr  i  te  t  5  /  97  ) 

Wr i t  e ( 1 , 96 ) 

Wr i te ( 1 / 97 ) 

Do  10  Index  =  1,4 

Urite(5,98)  QuadNo, XArray(l,  Index), YArray(l,  Index) 
Write(l,98>  QuadNo,  XArray(l,  Index), YArrayil,  Index) 
Cont  inue 
Wr i te ( 1 , 97 ) 

FoT-matiTS,  'Quadr  i  lateral ',  TIS,  '  X-Coor  d  ' ,  T30, 

' Y-Coord  '■  ) 

Forma  t(T3,  ' - ',  T18,  ' - ',  T30, 


90  Forma  t<  T9,  II,  T 1 3,  99  4,  T27,  F9.  4) 

Return 

End 


Subroutine  Trans formilnXcoord, InYcoord,  XArray, YArray, 

Ne'jiXArray,  NeuiYArray  ) 

This  subroutine  normaliics  the  inputted  qua d r  1 1  a t er a  1 3  to 
coordiriates  (0,0)  for  easy  calculations  based  on  the  origin. 

Real*8  Xarray(l,  4!',  VArray(l,4) 

RcalaQ  NeijjX Array!  i.,  4),  iNeuiYArrayf  1-  4) 

Rea  la-3  [nXcoord,  InYcoord 
Integer  Index 


Do  10  In OCX  =  1,4 

NeiiiXArrayi'l,  Index) 
MewYAr  -  a  y  (  1 ,  I  n  •1  a  x  ) 
Cont  inuc 


X Array! I,  Index) 
V  A  r  r  a  y ( 1 ,  Index) 

C-23 


I  n  X  c  o  o  r  t j 

InYcoord 


Return 

End 


•a-j:-* •#••«■■«■■)«•■)»•-»••<• -k-*-!*-** -is- •*■»•■«■  ■'(■■is- 
Subroutine  T e s t  _I  t  <  Ne'tiXAr ra y ,  NewYAr ra y  -  X >  Y  ) 

Inteqnr  EValuei 

R  eal  NeniXAr  r  ay  (  1  >  ^  .  NeiuYArray  (  1 .  4 ) 

Real ■'tH  X.  Y,  A(4)j  M<4).  Sum,  R,  Tmpl,  Tmp2 
Chara>:  ter-s-l  Loop, Exit 

R  =  0,  OoOl 

Do  10  Index  =  1,4 

M  (  i  nd  e  X  )  =  DSQR  T  i  (  Ne'xjX  Array  (1,  Index)  **  2)  + 

■t  ( NeujY  Array  (1,  Index)  ■»■«•  2)) 

If  <M(index)  .  LE.  R)  Then 
C  Pi-evention  of  division  by  zero. 

Wr  i  t  e  <  5,  ?S ) 

Wn  be  ('  1,  9S) 

GOTu  13 
Eiti  If 

li')  Conti-iiue 


J  =  2 
Sum  s  0 

Do  Wh  »  1  e  <  Ind  .?  X  .  I.E, 

If  V  I  nd  e  X  .  EO.  4  ) 

J  ^  1 
End  1  f 

finpl  =  ((NemXArrey 
( Ne(i;XAmy 
Tmn2  =  (  <  NeviiXArr ny 
( Ne'oYAr  ray  < 
A<  index)  =  DACOSui 
If  ( A ( I n  d  e  X )  .  G  T 
Wr  i  te  (  5<  30  ) 

Pormat(T2,  'The 
Wr  i  t  e  <  1 ,  3 1  ) 

Eor  ma  t  i  r2,  '  TIi  e 
GO  TO  13 
En  d  i  f 

If  (  i'mp  1  .  LT,  0  0 ) 


')) 

Ch  en 


(1.  Index)  *  Neti'YArr  ay  (  1 ,  J  )  ' 
i  1 ,  J )  *  NewY Array ( 1 ,  I nd  e  x ) )  ) 
<1,  Index)  ♦  Nen'XArray  (  1 ,  J  )  ) 
t ,  Index  )  *  Neu)Y Array  (  1 .  J  )  )  ) 
rmp2  /  (M( Index )  *  M( J) ) ) 
179, 88)  Then 

point  is  ON  THE  LINE.  ') 

point  is  ON  THE  LINE  ') 


.4  (  I  n>t  e  X  )  =  -A  (  T  ni 
End  1  f 

Sun)  ~  Sum  A'llnde 
I  O')  ex  -  Index  I 

J  =  +  1 

End  Du 

EV.ilue  =  JIDNNT(,>3 


Ti)  en 
lex  ) 


(  EV  J  1 !)  e  .  EQ. 
Wr  1 1  e  (  5.  27  ) 
Write( 1, 27) 


9 ( Sum  /  360  0 ) ) 
Ti)  en 


Wr 1 t  e ; 5 .  2S  > 
W  (•  1 1  e  I  1  •  29  ) 
End  i  f 


•  «  •  *.  ■  4  • 


V.  A  •'.  '  -  •  •  ^  - 

■  g.v  ^ 


27 

Forma  t ; 72- 

'Point 

not  enclosed.  ') 

20 

Format  i  T2.- 

'Point 

13  enclosed.  ') 

90 

Format ( T2- 

'Point 

of  interest  is  on  the  VERTEX') 

13 

Contiriuc 

R  e  t  u  T-  n 

End 

'oubrautine  SkipfLines) 

C 

C  This  subroutine  skips  the  indicating  number  of  lines 

C  the  place  inhere  skip  has  been  called  from  (i.e.  many 

/-» 

U 

Integer  I. Lines 

Do  10  I  =  1.  Lines 
Wr  i  ta  (  b,  11) 

10  Continue 

11  Forma t (  '  '  ) 

Return 
End 

C -4 -1- -a-* 4- 

C  •*-*■*•*-*■***•*■*  iS- ■*  T*  ^  -a--**  ■»•»•■»•»■»*•»■»-»■•* -s- •!*■■»■»■<•■<»••»•»■»-»  ■!»•■»■■»  ■» 

Subroutine  Promptl(Rasponse) 

C 

C  This  subroutine  Prompts  the  user  if  he/she  uants  to  input 

C  data  to  determine  if  the  point  of  interest  falls  within 

C  any  one  of  three  quadrilaterals. 

C 

Chara-:ter->l  Response 


by 

places) 


Call  Sk ip ( 1 ) 

Wr  i  te  (  5.  5.5  ) 

Wr  i  t  e  (  5<  5ij  ) 

Wr 1 te ( 5/  57 ) 

Call  Skip  a) 

Read(5/79)  Response 
11  Forma  t  (  ■'  '  ) 

55  Format(T2,  'Would  you  like  to  input  some  coordinates  to  '  ) 

5<3  Format(T2/  'determine  if  a  point  falls  in  an  inputted  ') 

57  FQrmar:iT2,  '  qua  d  r  i  1  a  t  e  r  a  1  Y/N  ') 

99  F o  r ma  t ( A  1  ) 

Return 

End 


Subroutine  Prompt2ine3ponse.  In X coord-  InYcoord) 

C 

C  This  su-b  routine  Anomots  the  user  if  he /  she  wants  to  input 

C  data  to  determine  if  the  point  of  interest  falls  within 

C  any  one  o-^  three  quadrilaterals,  c_25 


keal*8  InXcoord,  InYcoord 
Character*!  Response 


Call  SRip( I ) 
write(5,  55) 

Read(5,  89)  Response 
Call  8lcip(l) 

If  ((Response  .  EQ.  'Y')  .  OR.  (Response  .  EQ.  'y'))  Then 
write(  1  ,  01) 

Write(l,  01) 

WriCe( 5  ,  97  ) 

Read(3,  99)  InXcoord 
Write(l,  95)  InXcoord 
WriCe(l,  01) 

Wrice(5,  98) 

Read(3,  99)  InYcoord 
Write(l,  96)  InYcoord 
write(l,  01) 

Write(l,  01) 


01 

End  if 

F  0  rma  c ( ' 

') 

55 

Formac( T2 , 

'  Input  a  point  to  test??  Y/N  ') 

89 

torinac(  A1 ) 

95 

Formac(T2 , 

'Inputted  X-Coordinate  ===>  F9.6) 

96 

Format (T2 , 

'Inputted  Y-Coordinate  ===>  ',  F9.6) 

-il 

Format(T2 , 

'Input  X-coord:  ') 

98 

Format(T2 , 

'Input  Y-coord :  ' ) 

99 

Format( F9 . 

4) 

Return 

End 

C-26 


QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


C-1  0.0000  0.0000 

2.0000  0.0000 

A. 0000  0.0000 

6.0000  0.0000 


CHECKING  THE  INPUT  (X,  Y)  POINTS 


INPUT  X 


INPUT  Y  GENERATED  RESULT  ACTUAL  RESULT 


4.0000 

4.0000 

2.0000 

6.0000 

8.0000 

5.0000 

2.0000 

8.0000 

6.0000 

-3.0000 

-2.0000 

0.0000 

2.0000 

4.0000 

6.0000 

10.0000 

6.0000 

3.0000 

4.0000 

7.0000 

6.0000 

13.0000 

6.0000 

12.0000 

0.0000 

7.0000 


4.0000 

1.0000 

3.0000 

3.0000 

3.0000 

0.0000 

7.0000 

7.0000 

2.0000 

4.0000 

7.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

8.0000 

3.0000 

6.0000 

0,0000 

4.0000 

0.0000 

7.0000 

0.0000 

5.0000 

5.0000 


OUT 

OUT 

OUT 

OUT 

OUT 

ON  THE  LINE 
OUT 
OUT 
OUT 
OUT 
OUT 

CORNER 

CORNER 

CORNER 

CORNER 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 


OUT 

OUT 

OUT 

OUT 

OUT 

ON  THE  LINE 
OUT 
OUT 
OUT 
OUT 
OUT 

CORNER 

CORNER 

CORNER 

CORNER 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 

OUT 


QUADRILATERAL 

X-COORDINATE 

Y-COORDINATE 

C-2 

0.0000 

0.0000 

10.0000 

0.0000 

6.0000 

8.0000 

3.0000 

3.0000 

CHECKING  THE  INPUT  (X 

,  Y)  POINTS 

INPUT  X 

INPUT  Y 

GENERATED  RESULT 

ACTUAL  RESULT 

4.0000 

4.0000 

IN 

IN 

4.0000 

1.0000 

IN 

IN 

2.0000 

3.0000 

OUT 

OUT 

6.0000 

3.0000 

IN 

IN 

8.0000 

3.0000 

IN 

IN 

5.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

2.0000 

7.0000 

OUT 

OUT 

8.0000 

7.0000 

OUT 

OUT 

6.0000 

2.0000 

IN 

IN 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

0 .0000 

ON  THE  LINE 

ON  THE  LINE 

4.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

10.0000 

0.0000 

CORNER 

CORNER 

6.0000 

8.0000 

CORNER 

CORNER 

3.0000 

3.0000 

CORNER 

CORNER 

4.0000 

6.0000 

OUT 

OUT 

7.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

4.0000 

IN 

IN 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

IN 

IN 

12.0000 

0.0000 

OUT 

OUT 

0.0000 

5.0000 

OUT 

OUT 

7.0000 

5.0000 

IN 

IN 

QUADRILATERAL 

X-COORDINATE 

Y-COORDINATE 

C-3 

0.0000 

0.0000 

4.0000 

6.0000 

7.0000 

0.0000 

10.0000 

0.0000 

CHECKING  THE  INPUT  (X,  Y)  POINTS 


INPUT  X  INPUT  Y  GENERATED  RESULT  ACTUAL  RESULT 


4.0000 

4.0000 

IN 

IN 

4.0000 

1.0000 

IN 

IN 

2.0000 

3.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

3.0000 

OUT 

OUT 

8.0000 

3.0000 

OUT 

OUT 

5.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

2.0000 

7.0000 

OUT 

OUT 

8.0000 

7.0000 

OUT 

OUT 

6.0000 

2.0000 

ON  THE  LINE 

ON  THE  LINE 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

4.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

10.0000 

0.0000 

CORNER 

CORNER 

6.0000 

8.0000 

OUT 

OUT 

3.0000 

3.0000 

IN 

IN 

4.0000 

6,0000 

CORNER 

CORNER 

7.0000 

0.0000 

CORNER 

CORNER 

6.0000 

4.0000 

OUT 

OUT 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

OUT 

OUT 

12.0000 

0.0000 

OUT 

OUT 

0.0000 

5.0000 

OUT 

OUT 

7.0000 

5.0000 

OUT 

OUT 

QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


C-4 

0.0000 

0.0000 

6.0000 

4.0000 

13.0000 

0.0000 

6.0000 

0.0000 

CHECKING  THE  INPUT  (X,  Y)  POINTS 

INPUT  X 


INPUT  Y 


GENERATED  RESULT 


ACTUAL  RESULT 


4.0000 

4.0000 

2.0000 

6.0000 

8.0000 

5.0000 

2.0000 

8.0000 

6.0000 

-3.0000 

-2.0000 

0.0000 

2.0000 

4.0000 

6.0000 

10.0000 

6.0000 

3.0000 

4.0000 

7.0000 

6.0000 

13.0000 

6.0000 

12.0000 

0.0000 

7.0000 


4.0000 

1.0000 

3.0000 

3.0000 

3.0000 

0.0000 

7.0000 

7.0000 

2.0000 

4.0000 

7.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

8.0000 

3.0000 

6.0000 

0.0000 

4.0000 

0.0000 

7.0000 

0.0000 

5.0000 

5.0000 


OUT 

IN 

OUT 

IN 

OUT 

ON  THE  LINE 
OUT 
OUT 
IN 
OUT 
OUT 

CORNER 
ON  THE  LINE 
ON  THE  LINE 
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Integer  Guadcntr, I nd e i , Flag . NoOf 3 i d es 
Rcal»G  Sigma.  Temp.  X,  Y,  U(9.'#  V(8).  NoOfEnclosures 
Realms  S(3).  InX(8.P>.  I  riY  ( 3#  8 ) .  Tmp  1 8  >  #  De  1  ta  (  8 ) .  Ep  s  i  1  on 
Cnaractar^l  Answ. Resp 

Data  Quad  c  n  tr ,  Ep  3  i  1  on#  NoQ  PS  i  d  es  /O.O  0001.4/ 

Open(Unit  =  l.File='' Freds.  Dat'  #Statu3='New') 

Op  en  ( Un  i  t=2.  File='Pt5.Dat'  .Status='01d') 

0pen(Unit=4. File='FredRes.  Dat',Status='01d ') 

Op  en  (  Un  i  t=3.  File='Test.  Dat'  .Status='01d') 

Quadcntr  =  Guadcntr  +  1 
Index  =1 

Do  While  (Index  .  LE  NoOfSides) 

Wr i te ( 1 . 01 > 

W7-ite(5.99^  I  n  d  0  X  .  Guad  c  n  tr 
Wt  i  t  .a  (  1 .  99  )  I  n  d  e  .X .  Guadcntr 
Read (2. 97)  InX (Quadcntr# Index ) 

WTite(1.97)  InX(Quadcntr.  Index) 

Wr i ce ( 1  01  ) 

Writc(5#98)  Ind e X  .  Quad cntr 
Write(1.9a)  I  nd  e  .t ,  Quad  cn  tr 
Read(2. 97)  I nY ( Quad c n tr .  Ind e x ) 

Wr  i  t e ( 1 . 77 )  In Y (Quadcntr#  Index) 

Index  =  Index  +  1 
End  do 
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Wr i te (  1 ,  96 ) 

Read ( 0. 97 )  X 
Write  (1,97)  X 
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U(Quadcntr)  =  I n X ( Quad cn tr #  Index)  -  X 
Y'.Quadcntr)  =  I  nY  ( Quad  cn  tr ,  Index)  -  Y 
Sigma  =0 

Tmp(Quadcntr)  =  Scalc(U(Quadcntr).V^Quadcntr)) 

Flag  =  1 

I  n  d  e-  X  =2 

Wr  i  te  (  3,  31  )  X,  Y 

Formati72.  'Testing  point  C  ',F5  2.  '  '.F5  2.  ■  ] 

Do  While  ((Index  LE.  5)  AND  (Flag  EQ,  1)) 

If  (index  . EQ.  5 )  Then 
Index  =  1 
Flag  =  0 
End  1  f 

Ufl.ndex)  =  InX(Quadcntr#  Index)  -  X 
V (In. 3 ex)  =  In V (Quad cntr.  Index)  -  Y 

If  ((|J(  Index;  EQ.  0)  AND.  ('/(Index)  .  EQ.  O)  Tnon 


1  n  X  ( (juad  cn  tr  #  Index)  -  X 
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Scale  (U(Quadcntr).  '/(Quadcntr)  ) 
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Formati’73/  'Value 
Format<T2-  'Input 
For(nat(T3,  'Input 
FormatCF15.  4> 
Format(T2i  'Input 
Forma t(T2-  'Input 
Stop 
End 
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ON  THE  LINE 

6.0000 

4.0000 

IN 

IN 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

IN 

IN 

12.0000 

0.0000 

OUT 

OUT 

0.0000 

5.0000 

OUT 

OUT 

7.0000 

5.0000 

IN 

IN 

v;,Vj.Ny 


PSuGPAfI  HANCSZ'NESS  (  INPUT,  OUTPUT.  PT3.  TESiPTS.  MESS.  RESPONSES) , 

LAHEL  S9,  70i 

TYPE 


SETTINGS  =  RECORD 

ZERO 

INTEGER 

POS 

INTEGER 

NEO 

INTEGER 

NA 

INTEGER 

END; 

POINT  =  RECORD 

X,  Y 

REAL; 

FLAG 

PNT ARRAY 
HF ARRAY 
FIG 


END; 

(NEG,  POS,  ZERO,  NOTAPP); 

ARRAVCi.  43  OF  POINT; 
array: 1. .  43  OF  FLAG; 

(CONVEX,  CONCAUE,  TRIANG,  TRIRAY,  BOWTIE,  INVALID, 


LINE)  , 


VAR 

I,  J,  K 

CAGECONDITION 
WRONGS  I GN, QC, PC 
IMvEX 

POrNTi, PO INTO 
DETERM INANTVALUE 
MORI- WORK 
RE3P, AN3U 
3ETCHAR3 
SETCHARGOLD 
APGINT 
P 

HFLAG 

FIGURE, TESTFIG 
PTS, TESTPTS. MESS 
RESPONSES 


INTEGER; 

INTEGER; 

INTEGER; 

INTEGER; 

INTEGER. 

REAL; 

BOOLEAN; 

CHAR ; 

SETTINGS; 

SETTINGS, 

POINT; 

PNTARRAY; 

HFARRAY; 

FIG; 

TEXT; 

TEXT; 


PROCEDURE  DUMP  IHF  :  HFARRAY;  PNT  ;  PNTARRAY;  SETCHARS  :  SETTINGS.  PT  POINT; 
VAR 

I  :  INTEGER; 

BEGIN 

WR I TELN ('IN  DUMP  '  ) ; 


WRITELN(  ' - : - '  >  , 

WRITELN('  POINT  4  X-COORD  Y-COORD  HFLAG'), 

WRITELN(  ' - ' ) , 


FOR  I  :=  1  TO  4  DO 


4 


BEGIN 

WRITE('  ',1:2.'  ' .  P  C  I  3 ,  X .  6 .  3,  '  '.Pin  Y  6:3. 

WRITELN('  '.HFLAGCI3;, 

END; 

WRITELN(  ' - -  ) , 


UR  I  TELN  (' INPUT  X-COGRD  ==:>  PT  X.6  3,  '  INPUT  Y-COORD  =  =  :•  '  .  PT  Y  o  3). 

WR  I  TELN  ' -  I  . 

==>  '.SETCHARS  ZERO.  2), 

==>  ',  SETCHARS.  NEG:  2)  , 

==:•  '.SETCHARS  POS  2), 

==:>  SETCHARS  NA  2)  . 


PNTARRAY)  . 


C-44 


WRITELNt  'ZERO  SETTING 

WRITELNC  'NEGATIVE  SETTING 
UP ITELN-  'POSITIVE  SETTING 
WRITELN('NaT  APP  SETTING 

END, 

PROCEDURE  INITIAL! ZE^ VAR  PNT 
VAR 

I  :  INTEGER, 

BEGIN 


FUR  I  ; =  1  TD  4  DO 
BEGIN 

PNTlI3.  X  ;  =  0,  0; 
PNTc  13.  Y  :  =  0.  0; 

END; 

END; 


PROCEDURE  INir(UAR  HF  :  HFARRAY;  VAR  SETCHARS  :  SETTINGS); 
VAR 

I  :  INTEGER; 

BEGIN 

fur  I  :=  1  TO  4  DO 


HFCn  :=  NOTAPP; 


SETCHARS.  ZERO 

:  =  0, 

SETCHARS  POS 

:  =  0. 

SETCHARS.  NEG 

:  =  0; 

SETCHARS.  NA 

:  =  0, 

END; 


FUNCTION  SPECIFYFIG(HF  :  HFARRAY;  SETCHARS  ;  SETTINGS;  VAR  J  •  INTEGER)  F 
VAR 

I  :  INTEGER; 

BEGIN 

IF  < (SETCHARS.  NEG  =  4 )  OR  < SETCHARS.  POS  =  4))  THEN 
BEGIN 

.  SPECIFYFIG  :=  CONVEX 
END 
ELSE 

IF  < (SETCHARS.  ZERO  =  0)  AND 

( (SETCHARS.  POS  >  2)  OR  ( SETCHARS.  NEG  >  2)))  THEN 
BEGIN 

IF  (SETCHARS. NEG  =  3)  THEN 
BEGIN 

FOR  I  :=  1  TO  4  DO 

IF  HFCI3  =  POS  THEN 
J  :  =  I; 

END 

ELSE 

IF  (SETCHARS  POS  =  3)  THEN 
BEGIN 

FOR  I  .=  1  TO  4  DO 

IF  HFCn  =  NEG  THEN 
J  :  =  I. 

END 

ELSE 


WR I TELN ( ' DUMMY ' ) , 

SPECIFYFIG  CONCAVE, 

END 

ELSE 

IF  (SETCHARS  ZERO  =  1)  THEN 
3  EG  I N 

FOR  I  :=  1  TO  4  DO 

IF  HFCI3  =  ZERO  THEN 
J  :  ■=  I , 

IF  ((SETCHARS  NEG  =  2)  OR  ( SETCHARS.  PCS  =  2))  THEN 
SPECIFYFIG  :=  TRIRAY 
ELSE 

IF  ( (SETCHARS  NEG  =  3)  OR  ( SETCHARS.  PCS  =  3)'  THEN 
SPECIFYFIG  •=  TRIANG, 

END 


ELSE 
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Ir  c  (StTCHARS.  NE'3  =  2)  OR  < SETCHARS.  POS  =  2))  THEN 
SPECIFYEIG  :=  BGWTIE 
ELSE 

IF  (SETCHARS.  zero  =  4)  THEN 
SPECIFVFIG  :=  LINE 
ELSE 

SPECIFVFIG  :=  INVALID; 


FUNCTION  DETERNINANTCPl,  P2,  P3  :  POINT) 
VAR 

TMPl ,  TMP2.- TMP3  ,  TMP4  ;  REAL; 

BEGIN 


REAL. 


TMPl 

TMP2 

TMP3 

TMP4 


=  Pi.  X  ■*  (P2.  Y  -  P3.  Y)  ; 

=  P2.  X  ■*  (PI.  Y  -  P3.  Y> ; 

=  P3.  *  (PI.  Y  -  P2.  Y) ; 

=  TMPl  -  TMP2  +  TMP3i 


DETERMINANT  . =  TMP4; 


PROCEDURE  LOADPGINTCVAR  TEST  ;  POINT); 

BEGIN 

UiRITELN(  'ENTERING  LOADPQINT'  ); 

READLN(  TE5TPT3.  TEST.  X  ); 

URITELNCMESS,  'ENTER  X-COORD  TO  BE  TESTED  ===>  TEST  X  S  3 
READLN(  TESTPTS,  TEST  V  ). 

WRITELN(ME3S,  'ENTER  Y-CGORD  TO  BE  TESTED  ===;■  TEST  Y'6  3), 

WRITELN(  'LEAVING  LOADPOIMT'  ); 

END, 

PROCEDURE  LOaDFIGUR£(VAR  PNTS  :  PNTARRAY  ), 

VAR 

I  .  INTEGER, 

BEGIN 

FDR  I  :«  1  TO  4  DO 
BEGIN 

R£ADLN(  PT3,  PNT3CI]  X  ); 

URITELNi;  MESS,  'ENTER  X-COORD  OF  POINT  #',12,  '  ===;-  ',  PNTSC I  ]  X  ^  " 
READLN(  PT3,  PNTSCII.Y  ); 

WRITELNi  MESS,  'ENTER  Y-COORD  OF  POINT  #',I  2,  '  ===>  ■,  PMTSCI]  Y  S  3 
END, 


PROCEDURE  LOADHFLAGARRAY (  I 
BEGIN 

I,'"  DUAL  C  0.0  THEN 
HF  C  I  ]  :  =  NEG 
ELSE 

IF  DUAL  -  0.  0  T'HEN 
HFCn  :=  ZERO 
ELSE 

IF  DV^L  ;>  0.  0  THEN 
HFC  n  =  PCS 
ELSE 

HFC  n  :  =  NOTAPP, 


PROCEDURE  CALC'JLATESETTINGS!  J  INTEGER 
BEGE.N 

CASE  HFCJ]  OF 

POS  SC.  POS  .  =  3C  POS  -t-  I, 


INTEGER;  DVAL;  REAL-  VAR  HF  :  HFARRAY  ), 


INTEGER,  HF  HFARRAY,  VAR  SC  SETTINGS), 


mmmm 


NE3 
ZERO 
NOT APR 
END; 


30.  NEG 
SC.  ZERO 
SC.  NA 


=  SC.  NEG  +  1; 

=  SC.  ZERO  +  1; 
=  SC .  NA  +•  1  i 


FUNCi  rON  POINTINCONVEX (SC  :  SETTINGS)  ;  INTEGER; 
'  VAR 

FI.  F2,  F3  ;  BOOLEAN; 

BEGIN 


Ft 

: =  FALSE, 

F'.-' 

; =  FALSE, 

P-'i 

; =  FALSE, 

IF 

( (SC  PuS  =  4) 
FI  :  =  TRUE; 

OR 

(SC 

NEG  =  4) ) 

THEN 

T  ^ 

1  r 

<  (  SC .  PijS  =  3  1 
F2  :  =  TRUE; 

OR 

(SC 

NEG  =  3 ) ) 

THEN 

IF 

( ( SC  POS  =  2 ) 
F3  :  =  TRUE, 

OR 

( SC. 

NEG  =  2) ) 

THEN 

CASE  SC.  ZERO  OF 
0 


0  :  IF  Fi  THEN 

POINTINCONVEX  :=  1 
ELSE 

POINTINCONVEX  ;=  4. 

1  ;  IF  F2  THEN 

POINTINCONVEX  :=  2 
ELSE 

POINTINCONVEX  ;=  4; 

2  :  IF  F3  THEN 

POINTINCONVEX  : =  3 
ELSE 

POINTINCONVEX  :=  4; 
OTHERWISE  POINTINCONVEX  :=  4; 


FUNCTION  POINriNOTHER(VAR  I  :  INTEGER;  SC  :  SETTINGS)  BOOLEAN; 
VAR 

TMP  :  BOOLEAN, 

BEGIN 

TMP  FALSE, 

IF  ((SC.POS  >  0)  AND  <30.  NEG  >  0))  THEN 
I  :  =  4 
ELSE 

IF  SC.  IEhQ  =  2  THEN 
I  ;  =  T 
ELSE 

TMP  -  TRUE, 

Pi  (INTINOTHER  =  TMP; 

END; 

PROCEDURE  COMPARE; A,  3  SETTINGS;  VAR  J  :  INTEGER), 

BEGIN 

ir  (((A.  ZERO  •  0)  AND  ((A  POS  =  3)  OR  ( A.  NEG  =  3))) 

OR  ((A.  ZERO  -  1)  AND  (  li  ZERO  =  1)))  THEN 

INDEX  =  4 
CL3E 

IF  A.  ZERO  ;  2  THEN 
INDEX  J 
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ir  i(A.  ZERO  3.  ZERO  =  i)  AND 

i(A.  PQS  =  2)  OR  (A.  NEG  =  2)))  THEN 
INDEX  : =  2 
ELSE 

INDEX  : =  1; 

EMDi 

FUNCTION  POINTINTRlANGScSC  :  GETTINGS)  :  INTEGER; 

B  E'^  I  tv 

ir  ((SC.  ZERO  =  1)  AND  ((SC.  PQS  =  2)  OR  (SC.  NEG  =  2)))  THEN 
POINTINTRIANGS  :=  2 
t-LSE 

Ir  (SC.  ZERO  =  2)  THEN 
POINTINTRIANGS  ;=  3 
ELSE 

IF  (C'SC.POS  -  3)  OR  (SC.  NEG  =  3))  AND  (SC.  ZERO  =  0?)  THEN 
POINTINTRIANGS  :=  1 
ELSE 

POINTINTRIANGS  : =  4; 

END; 

PROCrOURE  OUTPUTNESSAGEC INDEX  :  INTEGER;  P  :  POINT); 

BEGIN 

t'ASE  INDEX  OF 

1  :  BEGIN 

WRITELNi  'THE  POINT  C P,  X:  6:  3,  .  P.  Y :  6:  3,  • '  ]  IS  INTHEFIGuRE 

WRITE  (MESS,  'THE  PO  I  NT  C  ' .  P .  X :  6 ;  3.  ' ,  P.  Y:  6;  3)  ; 

WRITELN(MESS,  ']  IS  IN  THE  FIGURE.  '  ) ; 

END; 

2  :  BEGIN 

WRITE(  'THE  POINT  C  P.  X;  6:  3,  ',P.Y:6:3,  ']  ISON'), 

WRITELN('  THE  BOUNDARY  OF  THE  FIGURE.  '); 

WRITE(MESS,  'THE  POINT  C P.  X:  6. 3,  ' ,  P .  Y  •  6 .  3,  '  3  ISON), 

WRITELNCMEES,  '  THE  BOUNDARY  OF  THE  FIGURE.  '); 

END; 

3  :  BEGIN 

WRrTE(  'THE  POINT  C ',  P.  X;  6:  3,  ',  ',P.Y:6:3,  '3  ISON'), 

WfilTELNf'  THE  VERTEX  OF  THE  FIGURE.  ■' ) ; 

WRITE(ME33,  'THE  POINT  C P  X:  6:  3,  ,  P.  Y:  6:  3,  '  3  ISON'), 

WRITELN(MEG3,  '  THE  VERTEX  OF  THE  FIGURE.  '), 

END, 

4  :  BEGIN 

WRITE(  'THE  POINI  C  ' ,  P.  X;  6;  3,  ' ,  ',P,Y:6:3,  '3  IS  OUT  '  ) , 

WRITELN('Or  THIS  FIGURE.  '), 

WRITE(MEGS,  'THE  POINT  C',P  X:6  3,  ',  ',P.Y:6;3,  '3  IS  OUT  '), 

WRITELN(M£SS,  'OF  THIS  FIGURE.  '); 

END; 

OTHERWISE  WR ITELNCMESG,  'ERROR'), 

END; 

END; 

P.POC,~DURE  CALC  DETER."!  INANT  i  VA.R  H.'^  ,  HFARRAY,  VAR  SC  :  SETTINGS,  I  INTEGER, 

PI,  P2,  P3  :  POINT,  VAR  DVAL  REAL  ), 

DEG  O'! 

L>'.''AL  =  DETL.RM  inant  .  P  ;  ,  ;-2,  P3  ) ; 

LOAnHFLAGA.Rr'.A'*'  (  I,  OVAL,  HD  , 

C  ALt-  ULATElE  '  i  1  NG^  t.  I ,  HE ,  lC  .'  , 

END, 

BEGIN  ( MAIiN 

ursp  =  'V-,  C_48 


I TE  (  MK33 ) ; 

R!-SET<PTS)i 

Kt-:SET(RE:3P0NS(-:3;i 

REPEAT 

QC  :  ■  ’iC  +  1  j 
RESET (TEPTPTS)- 
PC  :  0; 

INir(AI.IZE<P)i 
I.OADKIGURE\P  ); 

REPEAT 

PC  :  ^  PC  -t-  1; 

WHI  i'ELN(  ' - '  )  . 

WHi  rELN(MEE;3) , 

UhlTELN'iMESS.  'CHECKING  QUADRILATERAL  #  QC:2  ), 

WKl iELN(  'CHECKING  QUADRILATERAL  #  ',  QC ;  2  )> 

WHI  rELN<  'CHECKING  POINT  «  PC:2  ); 

WK I  1 ELN  V  ' - '  ) , 

INIT(HFLAC. 3E I  CHAPS); 

I.NADPQINT;  APUINT  ); 

Of-  OUMP(HFLAG,  P.  3ETCHAR3.  APOINT),  *) 

FOR  I  :  =  1  rO  4  DO 
BEGIN 

J  :  =  C  +  1. 

IF  J  -  5  THEN 
J  :  =  1; 

K  :  ■=  J  +  li 
IE  K  -  3  THEN 
K  ;  =  li 

CALCD£TERHINANT<HFLAG,  SETCHARS,  I,  PCn>  P[JD,  PCKD, 

DETERM I NANTV ALUE ) , 

END.- 

(■>DUMH<HF1.AG,  H,  SE  i  CHARS.  APOINT);  *) 

WRONGS I QN  :  -  55; 

TE3TFIG  :  ^  SPEC  I FYF IG < HFLAG. SETCHARS, WRONGSI GN ) , 

WRONGS I GN  ;  -  WRONGS I GN  +  1; 

POINT  I  :  =  WRONGS I ON; 

IE  WRONGS  I GN  >=  3  THEN 
P0INT2  : =  WRONGS I GN  -  2 
El  SE 

P0INT2  : =  WRONGSIQN  +  2; 

CASE  TESTFIG  OF 
CONVEX  BEGIN 

WRIfELNI  'CONVEX  CASE'  ), 

Fi)R  I  :=  1  TO  4  DO 
BEGIN 

J  ;  =  I  I, 

IF  U  =  5  THEN 
J  :  =  1, 

CALCDETERMINANTCHFLAG,  SETCHARS.  I .  P  C  I  ]  .  .=  £  J 
APOINT,  DETERMINANT value ’ . 


END; 

OU  iPUTMESSAGEIPOINTlNCONVEX  (  SETCHARS  /  .  APOINT  , 


END; 

CONCAVE  BEGIN 


WR I TELNI 'CONCAVE  CASE'), 

IE  WRONCSIGN  =  5  THEN 
WRONGS  I  GN  .=  1, 

OUMH I HFLAG, P, SETCHARS. APOINT ) , 
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TRIRAY, 
TR lANG 


IMI I (HFLAG/ bEiuHARb); 

FOR  I  ;=  1  TO  3  DO 
BEGIN 

WRONGSIGN  :=  WRQNG3IGN  +  1; 

IF  WRONGSIGN  =  5  THEN 
WRONGSIGN  :=  1; 

K  ;=  (WRONGSIGN  MOD  4)  +  1; 

IF  I  =  3  THEN 

K  ; =  ( K  MOD  4 )  +  1; 

(*DUMP(HFLAG,  P,  3ETCHAR3,  APGINT  /  , 
CALCDETERMINANTtHFLAG,  SETCHAR3-  I, 

PCWRONG3IGN3, PlK], APGINT, 
DETERM I NANTVALUE  > , 

END; 

DUMP(HFLAG. P- 3ETCHARS, APOINT) ; 

MORFWORK  ;=  POINTINOTHER ( INDEX, SETCHAR5 ) , 

WR  ITEi  'BEFORE  MQROWORK.  MOREWORK  ==>  ' ,  MOREWLR':*. ; 
WRITELN(  '  INDEX  ==>  INDEX:  2  ), 

IF  NOT  MOREWORK  THEN 

QUTPUTME3SAGE( INDEX,  APOINT) 

ELSE 

BEGIN 

WRONGSIGN  :=  POINT!; 

IF  WRONGSIGN  =  5  THEN 
WRONGSIGN  :=  1; 

SETCHARSOLD  : =  SETCHARS, 

DUMP(HFLAG,  P,  SETCHARS, APOINT) , 

INITCHFLAG. SETCHARS); 

K  :=  WRONGSIGN; 

FOR  I  :=  1  TO  3  DO 
BEGIN 

IF  I  =  1  THEN 

WRONGSIGN  :=  WRONGSIGfJ  -  1, 

IF  WRONGSIGN  =  0  THEM 
WRONGSIGN  : =  4, 

K  :=  (WRONGSIGN  MOD  4)  1, 

IF  I  =  3  THEN 

K  ;  =  ( K  MOD  4 )  i , 

(*DUMP(HFLAG,  P,  SETCHARS,  aPQIMT)  .  ) 

CALCDETERMINANT(HFLAG, SETCHARS.  I. 

P  C  WRONGS  IGN  ]  .  P  CK  ]  ,  APO  I Ni" , 
DETERM I NANTVALUE  > , 
WRONGSIGN  ;=  WRONGSIGN  +  1, 

IF  WRONGSIGN  =  5  THEN 
WRONGSIGN  ;=  1, 

END; 

DUMPiHFLAG,  P.  SETCHARS, aPOINT) , 

MOREWORK  =  POINTINaTHER( INDEX, SETCHARS) . 
COMPAREiSETCHARS, SETCHARSOLD, INDEX), 
QUTPUTMESSAGE( INDEX, APOINT) , 


BEGIN 

WRITELN( 'TRI  CASE’ >, 

INITCHFLAG, SETCHARS) , 

FOR  I  •=  I  TO  3  DO 
BEGIN 

WRONGSIGN  :=  WRONGSIGN 
IF  WRONGSIGN  =  5  THEN 
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WRONGS  I GN  ;=  1, 

K  ;=  (WPONGSIGN  MOD  4)  -*■  1; 

IF  I  =  3  THEN 

K  :  =  ( K.  MOD  4 )  +  1 ; 

<*DUMP (HFLAG,  P,  SETCHARS.  AROINT  > ;  *) 
CALCDETERMINANKHFLAG,  SETCHARS-  I- 

PCWRONGSIGN], rCK], AFQ 
DETERM I NANTVALUE ) ; 

ENDi 

DUMP (HFLAG- P, SETCHARS- APOINT)- 

OUTP  UTMESS AGE (POINTINTRIANGS(SETCHARS)-APOINT 

ENDj 

INVALID, 

QCWT IE- 

LINE  ■  BEGIN 

WR I rELN( 'INVALID  CASE ' ) ; 

GOTO  69; 

El'-iD- 

END- 

READLN< RESPONSES- ANSW); 

WK I TELN'  'ANOTHER  POINT??  Y/N'.ANSW); 

GOTO  70, 

ANSW  =  'N  '  - 

WRITELN<  'THIS  FIGURE  “  ' , TESTFIG.  ' "  IS  AN  ILLEGAL  FIGURE 
WH ITELN(  'PLEASE  INPUT  ANOTHER  FIGURE  OR  TERMINATE  '). 

UNTII.  ((ANSW  =  'N')  OR  (ANSW  =  'n')); 

CLOSE (TESTPTS); 

READI  N  (  RESPONSES •  RESP  ) ; 

WR I TELN( 'ANOTHER  QUAD??  Y/N',  RESP); 

UNTIL  ((RESP  =  'N')  OR  (RESP  =  'n')); 

WRIT£LN(  'THAT''S  ALL  FOLKS!!'); 

CLOSE (PT3) i 
CLOSE (RESPONSES); 

CLOSE (MESS) - 
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QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


0.0000 

2.0000 

4.0000 

6.0000 
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LINE  IS  ILLEGAL  FIGURE,  NO  POINT  TESTED 


0.0000 

0.0000 

0.0000 

0.0000 


INPUT  X 


INPUT  Y 


GENERATED  RESULT 


ACTUAL  RESULT 


4.0000 

4.0000 

OUT 

4.0000 

1.0000 

OUT 

2.0000 

3.0000 

OUT 

6.0000 

3.0000 

OUT 

8.0000 

3.0000 

OUT 

5.0000 

0.0000 

ON  THE  LINE 

2.0000 

7.0000 

OUT 

8.0000 

7.0000 

OUT 

6.0000 

2.0000 

OUT 

-3.0000 

4.0000 

OUT 

-2.0000 

7.0000 

OUT 

0.0000 

0.0000 

CORNER 

2.0000 

0.0000 

CORNER 

4.0000 

0.0000 

CORNER 

6.0000 

0.0000 

CORNER 

10.0000 

0.0000 

OUT 

6.0000 

8.0000 

OUT 

3.0000 

3.0000 

OUT 

4.0000 

6.0000 

OUT 

7.0000 

0.0000 

OUT 

6.0000 

4.0000 

OUT 

13.0000 

0.0000 

OUT 

6.0000 

7.0000 

OUT 

12.0000 

0.0000 

OUT 

0.0000 

5.0000 

OUT 

7.0000 

5.0000 

OUT 
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QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


i 


r.i 
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I  t‘-. 
v!'. 


i 


I 

K 

0 


& 


*  . 


INPUT  X 


4.0000 

4.0000 

2.0000 

6.0000 

8.0000 

5.0000 

2.0000 

8.0000 

6.0000 

-3.0000 

-2.0000 

0.0000 

2.0000 

4.0000 

6.0000 

10.0000 

6.0000 

3.0000 

4.0000 

7.0000 

6.0000 

13.0000 

6.0000 

12.0000 

0,0000 

7.0000 


INPUT  Y 

GENERATED  RESULT 

ACTUAL  RESULT 

4.0000 

IN 

IN 

1.0000 

IN 

IN 

3.0000 

OUT 

OUT 

3.0000 

IN 

IN 

3.0000 

IN 

IN 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

7.0000 

OUT 

OUT 

7.0000 

OUT 

OUT 

2.0000 

IN 

IN 

4.0000 

OUT 

OUT 

7.0000 

OUT 

OUT 

0.0000 

CORNER 

CORNER 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

0.0000 

CORNER 

CORNER 

8.0000 

CORNER 

CORNER 

3.0000 

CORNER 

CORNER 

6.0000 

OUT 

OUT 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

4.0000 

IN 

in" 

0.0000 

OUT 

OUT 

7.0000 

IN 

IN 

0.0000 

OUT 

OUT 

5.0000 

OUT 

OUT 

5.0000 

IN 

IN 

QUADRILATERAL 


X-COORDINATE 


y-COORDINATE 


C-3 

0.0000 

0.0000 

4.0000 

6.0000 

7.0000 

0.0000 

LO.OOOO 

0.0000 

CHECKING  THE  INPUT  (X,  Y)  POINTS 

INPUT  X  INPUT  y 

GENERATED  RESULT 

ACTUAL  RESULT 

4.0000 

4.0000 

IN 

IN 

4.0000 

1.0000 

IN 

IN 

2.0000 

3.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

3.0000 

OUT 

OUT 

8.0000 

3.0000 

OUT 

OUT 

5.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

2.0000 

7.0000 

OUT 

OUT 

8.0000 

7.0000 

OUT 

OUT 

6.0000 

2,0000 

ON  THE  LINE 

ON  THE  LINE 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

4.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

10.0000 

0.0000 

OUT* 

CORNER 

6.0000 

8.0000 

OUT 

OUT 

3.0000 

3,0000 

IN 

IN 

4.0000 

6.0000 

CORNER 

CORNER 

7.0000 

0.0000 

CORNER 

CORNER 

6.0000 

4.0000 

OUT 

OUT 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

OUT 

OUT 

12.0000 

0.0000 

OUT 

OUT 

0.0000 

5.0000 

OUT 

OUT 

7.0000 

5,0000 

OUT 

OUT 

♦THIS  CORNER  IS  DEGENERATED  TO  ZERO  ANGLE.  THIS  ANSWER  DEPENDS  ON  DEFINITION. 


QUADRILATERAL 


X-COORDINATE 


y-COORDINATE 


C-4  0.0000  0.0000 

6.0000  4.0000 
13.0000  0.0000 
6.0000  0.0000 


CHECKING  THE  INPUT  (X,  Y)  POINTS 


INPUT  X  INPUT  Y  GENERATED  RESULT  ACTUAL  RESULT 


4.0000 

4.0000 

OUT 

OUT 

4.0000 

1.0000 

IN 

IN 

2.0000 

3.0000 

OUT 

OUT 

6.0000 

3.0000 

IN 

IN 

8.0000 

3,0000 

OUT 

OUT 

5.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

2,0000 

7.0000 

OUT 

OUT 

8.0000 

7,0000 

OUT 

OUT 

6.0000 

2.0000 

IN 

IN 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

4.0000 

0,0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

0.0000 

CN  THE  LINE* 

CORNER 

10.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

8,0000 

OUT 

OUT 

3.0000 

3.0000 

OUT 

OUT 

4.0000 

6,0000 

OUT 

OUT 

7.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

4.0000 

CORNER 

CORNER 

13.0000 

0.0000 

CORNER 

CORNER 

6.0000 

7.0000 

OUT 

OUT 

12.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

0.0000 

5.0000 

OUT 

OUT 

7.0000 

5.0000 

OUT 

OUT 

*THIS  CORNER  IS  DEGENERATE  AND  MAY  BE  CONSIDERED  A  LINE.  THIS  ANSWER 
DEPENDS  ON  DEFINITION. 


QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


C-5 

0.0000 

0.0000 

6.0000 

7 .0000 

12.0000 

0.0000 

6.0000 

3.000O 

CHECKING  THE  INPUT  (X,  Y)  POINTS 

INPUT  X 


iW5 


INPUT  Y 


GENERATED  RESULT 


ACTUAL  RESUL 


4.0000 

4.0000 

IN 

IN 

4. 0000 

1 .0000 

OUT 

OUT 

2.0000 

3.0000 

OUT 

OUT 

6.0000 

3.0000 

CORNER 

CORNER 

a. 0000 

3.0000 

IN 

IN 

5.0000 

0.0000 

OUT 

OUT 

2.0000 

7.0000 

OUT 

OUT' 

s.oooo 

7.0000 

OUT 

OUT 

6.0000 

2.0000 

OUT 

OU'I 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

0.0000 

OUT 

our 

4.0000 

0.0000 

OUT 

OUT 

6.0000 

0.0000 

OUT 

OUT 

10.0000 

0.0000 

OUT 

OUT 

6.0000 

8.0000 

OUT 

OUi' 

3.0000 

3.0000 

IN 

IN 

4.0000 

6.0000 

OUT 

OUT 

7.0000 

0.0000 

OUT 

OUT 

6.0000 

4.0000 

IN 

IN 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

CORNER 

CORNER 

12.O000 

0.0000 

CORNER 

Corner 

0.0000 

5.0000 

OUT 

OUT 

7 .OOOO 

5.0000 

IN 

IN 

QUADRILATERAL 

X-COORDINATE 

Y-COORDINATE 

C-6 

0.0000 

0.0000 

0.0000 

5.0000 

7.0000 

0.0000 

7.0000 

5.0000 

CHECKING  THE  INPUT  (X,  Y)  POINTS 
BOV«/TIE  IS  AN  ILLEGAL  FIGURE.  NO 

POINTS  TESTED 

INPUT  X 

INPUT  Y 

GENERATED  RESULT 

ACTUAL  RESULT 

4.0000 

4.0000 

OUT 

4.0000 

1.0000 

OU'L' 

2.0000 

3.0000 

IN 

6.0000 

3.0000 

IN 

8.0000 

3.0000 

OUT 

5.0000 

0.0000 

OUT 

2  .0000 

7.0000 

OUT 

8.0000 

7.0000 

OUT 

6.0000 

2.0000 

IN 

-3.0000 

4.0000 

OUT 

-2.0000 

7.0000 

OUT 

0.0000 

0.0000 

CORNER 

2.0000 

0.0000 

OUT 

4.0000 

0.0000 

OUT 

6.0000 

0.0000 

OUT 

10.0000 

0.0000 

OUT 

6.0000 

8.0000 

OUT 

3.0000 

3.0000 

OUT 

4.0000 

6.0000 

OUT 

7.0000 

0.0000 

CORNER 

6.0000 

4.0000 

IN 

13.0000 

0.0000 

OUT 

6.0000 

7.0000 

OUT 

12.0000 

0.0000 

our 

U.OOOO 

5.0000 

CORNER 

7.0000 

5.0000 

CORNER 

QUADRILATERAL 


X-COORDINATE 


Y-COORDINATE 


0.0000 

10.0000 

6.0000 

2.0000 


0.0000 

0.0000 

8.0000 

4.0000 


CHECKING  THE  INPUT  (X,  Y)  POINTS 


INPUT  X 


INPUT  Y 


GENERATED  RESULT 


ACTUAL  RESULT 


4.0000 

4.0000 

IN 

IN 

4.0000 

1.0000 

IN 

IN 

2.0000 

3.0000 

IN 

IN 

6.0000 

3.0000 

IN 

IN 

8.0000 

3.0000 

IN 

IN 

5.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

2.0000 

7.0000 

OUT 

OUT 

8.0000 

7.0000 

OUT 

OUT 

6.0000 

2.0000 

IN 

IN 

-3.0000 

4.0000 

OUT 

OUT 

-2.0000 

7.0000 

OUT 

OUT 

0.0000 

0.0000 

CORNER 

CORNER 

2.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

4.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

10.0000 

0.0000 

CORNER 

CORNER 

6.0000 

8.0000 

CORNER 

CORNER 

3.0000 

3.0000 

IN 

IN 

4.0000 

6.0000 

ON  THE  LINE 

ON  THE  LINE 

7.0000 

0.0000 

ON  THE  LINE 

ON  THE  LINE 

6.0000 

4.0000 

IN 

IN 

13.0000 

0.0000 

OUT 

OUT 

6.0000 

7.0000 

IN 

IN 

12,0000 

0.0000 

OUT 

OUT 

0.0000 

5.0000 

OUT 

OUT 

7,0000 

5.0000 

IN 

IN 

C-58 

.  '  V  *  W  •  L  *  W  »  ■  •  ^  ^ 

wV 


I 


V'.' 

5 


u 


L-'' 

It 


K.'. 

U 


|i 

I 

E 


s 


APPENDIX  D 

SOME  MATHEMATICAL  PROOFS 

Proposition:  The  point  P  lies  within  the  triangle  ABC  if,  and  only  if,  the 
extension  of  the  line  segments  AP,  BP,  and  CP  intersect  BC,  CA,  and  AB, 
respectively. 

Proof:  First  note  that  for  any  two  points  X  and  Y,  XY\X,  and  XY\Y  are 
connected  sets  (Moore,  Theorem  98).  Suppose  P  lies  within  ABC,  then  so  do 
AP,  BP,  and  CP  (by  Theorem  33,  Moore  applied  to  AP\A,  etc).  Since  AP 
divides  the  angle  CAB,  its  extension  intersects  BC,  and  similarly  for  BP 
and  CP. 

Suppose  one  of  the  line  extensions  does  not  intersect  the  appropriate 
side.  Without  loss  of  generality,  suppose  the  line  passing  through 
A  and  P  does  not  intersect  BC.  Since  it  does  not  pass  through  C,  PAC  are 
not  colinear,  and  similarly  PAB  are  not  colinear.  Thus  it  intersects  ABC 
at  the  endpoint  A  only.  Thus  the  connected  set  AP\A  does  not  intersect 
ABC,  and  hence  falls  entirely  outside  ABC  (again  by  Moore,  Theorem  33). 
Therefore  P  is  not  enclosed  by  ABC. 


y  / 


Rotation  of  coordinates: 

See  Figure  D-1  (a) 

X'  =  (x  +  y  tan  6)0036 
=  X  cos  6  +  y  sin 6 
y'  =  y/cos6  -  (x  +  y  tan  6)sine 

=  -X  sin6  +y(1/cos6  -  sir? S' cos  6) 
=  -X  sin  6  +  y  cos  6 

ic*  =  (-X  cose  +  y  sin  6) 
ir'  =  (-X  sine  +  y  cos  6) 


ANOTHER  PROOF 


Given;  three  ordered  points  Pi(xi,y^), 

yi  1 

det  *2  72  1  = 

*3  ^3  1 

Prove: 

h  =  0  if  P^  is  on  the  extended  line  P^  P2. 
h  >0  if  P3  is  on  the  left  side  of  the  directed  line  segment, 
h  <0  if  P^j  is  on  the  right  of  the  directed  line  segment. 

?r<i<3Z.i 

Translate  origin  to  P^  (see  Figure  D-1  (b)) 

Xi  =  (Xi  -x^)  where  (Xi,Yi)  =  (0,0). 

^i  »  (7i  -yi). 

Rotate  the  axis  so  that  the  X2  lies  on  the  axis,  (see  Figure  D-1  (c)) 
sine  =  ~'[2/3qrt(X2^*Y2^'^ 

COS0  =  X2/sqrt(X2^+Y2^)- 

Y3'  s  -X3  sine  ♦  Y3  cose 

Turn  in  transition  along  Pi,P2,p3,  is  left  if  Y’3  is  >0,  along  a  lii 


there  is  no  turn  if  Y’^  =  0,  and  turn  is  right  if  y'3<0.  If  this  function 
is  called  sign  then: 

1*3  s  -X3  ain6  +  Y3  cosG 

=  -X3X2/aqrt(X2^+Y2^^  X3X2/sqrt(X22+y22) 

s  [-X3Y2  +  X2Y3]/sqrt(X2^+Y2^) 

=  [-(X2-x-|)(  72-71 )  +  (x2-Xi)(73-7i)]/sqrt(X2^+X2^^ 
s  [  X371+X172+X273+X1 71-X372-X1 71  “X27i  -xi  73 1  /  sqrt  ( • 

yi  1 

sign  (det  X2  y2  1  )  =  sign  h  q.e.d. 

X3  73  1 

since  the  square  root  is  alwa7s  positive. 


